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Abstract 

The  Global  Positioning  System  (GPS)  constellation  provides  ranging  information 
that  delivers  inexpensive,  high  precision  positioning  for  terrestrial  users.  Satellites 
in  Low  Earth  Orbit  (LEO)  can  use  an  onboard  GPS  receiver  resulting  in  meter-level 
navigation  solution  accuracy.  There  are  limitations  to  using  GPS  for  positioning  for 
satellites  above  LEO.  The  number  of  GPS  satellites  whose  signal  can  be  received  de¬ 
creases  as  the  receiver’s  altitude  approaches  that  of  the  GPS  constellation.  Above  the 
GPS  constellation,  the  available  GPS  signals  for  ranging  will  originate  from  satellites 
on  the  opposite  side  of  Earth.  This  research  calculates  the  available  GPS  signals  to 
the  receiver  and  determines  the  expected  position  error,  while  considering  the  effects 
from  low  signal  to  noise  ratios,  poor  geometry,  and  signal  shift  caused  by  high  relative 
velocity. 
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EXPECTED  POSITION  ERROR  FOR  AN  ONBOARD  SATELLITE  GPS 


RECEIVER 

I.  Introduction 

1.1  GPS  Based  Orbit  Determination 

The  first  NAVSTAR  Global  Positioning  System  (GPS)  was  launched  in  1978  and 
has  quickly  become  the  leading  tool  for  Position,  Navigation,  and  Timing  (PNT) 
applications  worldwide.  GPS  receivers  have  become  relatively  inexpensive  and  ter¬ 
restrial  users  of  GPS  are  capable  of  obtaining  centimeter  level  accuracy.  The  low 
cost,  reliability,  and  success  rate  of  GPS  for  PNT  on  Earth  drove  the  desire  for  space 
based  application.  Benefits  of  implementing  the  already  operational  GPS  system  for 
on-board  orbit  determination  are  listed  below  [3]: 

1.  Reduced  amount  of  ground  stations  required  for  operation. 

2.  Scheduling  of  ground  station  operation  data  collection  made  easier  and  can  be 
done  in  advance 

3.  Reduced  cost  from  the  reduction  of  ground  stations 

4.  Orbit  determination  automation  possible 

5.  Multiple  satellite  formation  flying  possible 

When  the  receiver  remains  close  to  the  Earth’s  surface,  there  is  little  change  to 
the  expected  GPS  coverage  compared  to  a  ground  user.  GPS  signals  still  originate 
from  above,  Earth  obstructs  all  the  signals  from  underneath,  and  the  ionosphere  still 
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affects  the  signal.  Depending  on  the  altitude,  the  user  may  not  be  affected  by  the 
troposphere  but  there  are  other  consumer  based  challenges,  such  as  the  requirement 
for  an  unlocked  GPS  receiver  plus  the  issue  of  handling  increased  Doppler  magnitudes 
compared  to  ground  users.  Concerns  about  the  GPS  coverage  environment  increase 
when  the  receiver’s  altitude  approaches  or  surpasses  the  orbital  altitude  of  the  GPS 
constellation.  When  planning  for  GPS  utilization  at  or  above  the  GPS  orbital  altitude, 
issues  such  as  lower  signal  strength  from  increased  range,  error  sensitivity  due  to  poor 
geometry,  and  effects  of  increased  GPS  satellite-to-receiver  relative  velocity  have  to 
be  considered. 

1.2  Problem  Statement 

The  GPS  antenna  is  designed  with  Earth  coverage  as  the  main  focus  so  it  uses  a 
directional  beam  antenna  to  focus  its  transmit  power  over  Earth’s  hemisphere.  The 
area  of  this  ’main  beam’  is  located  between  GPS  and  Earth,  and  has  some  spillage 
over  Earth’s  edge.  The  availability  of  GPS  Space  Vehicles  (SVs)  decreases  as  a  user 
approaches  and  even  surpasses  the  GPS’s  constellation  altitude  due  to  being  outside 
of  the  GPS  main  beam’s  area.  Instead  of  relying  on  GPS  signals  that  lie  between 
the  GPS  SV  and  Earth,  users  above  the  GPS  constellation  depend  on  the  main 
beam  spillage  originating  from  the  opposite  side  of  Earth.  Multiple  issues  need  to 
be  reviewed  when  relying  solely  on  GPS  signals  originating  from  the  opposite  side 
of  Earth.  These  issues  include  limitations  on  the  number  of  available  satellites  due 
to  Earth  obstruction;  increased  error  caused  by  lower  signal  strength;  and  decreased 
GPS  angle  of  separation  relative  to  the  user,  leading  to  poor  geometry.  Conventional 
tools  for  position  determination  and  position  error  sensitivity,  such  as  the  Newton- 
Raphson  Least  Squares  algorithm  and  the  Dilution  of  Precision  (DOP)  metric,  are 
evaluated  for  applicability  in  these  scenarios. 
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1.3  Research  Contributions 


The  goal  of  this  research  effort  is  to  investigate  the  error  expected  when  calcu¬ 
lating  position  with  a  GPS  receiver  onboard  a  subject  SV  in  seven  different  orbital 
altitudes.  Position  error  is  estimated  at  each  point  in  the  desired  satellite’s  orbit. 
The  orbits  of  interest  for  calculating  position  error  include  those  located  below  as 
well  as  above  the  GPS  constellation,  for  example  Low  Earth  Orbit  (LEO)  and  Geo¬ 
stationary  Orbit  (GEO).  This  research  will  take  into  account  User  Equivalent  Range 
Error  (UERE),  signal  strength  degradation  due  to  distance,  and  the  effects  of  Doppler 
resulting  from  the  high  velocity  of  the  spacecraft  near  perigee.  Effects  of  the  GPS 
antenna  sidelobes  will  also  be  considered  when  determining  GPS  signals  available  to 
the  user  at  each  time  step.  This  thesis  outlines  a  process  for  evaluating  expected  GPS 
accuracy  tailored  to  the  concerns  of  satellite  orbits. 

1.4  Thesis  Overview 

The  purpose  of  this  research  is  to  evaluate  the  magnitude  of  errors  observed  when 
using  GPS  for  positioning  at  different  select  orbits.  To  better  understand  the  problem, 
a  look  at  orbit  characteristics,  the  GPS  constellation,  existing  errors  affecting  GPS, 
existing  position  determination  algorithms,  and  other  relevant  topics  are  presented 
in  Chapter  II.  Chapter  II  ends  with  the  introduction  of  the  resources  used  in  this 
research.  The  methodology  behind  the  construction  of  the  simulation  to  investigate 
this  problem  and  all  underlying  assumptions  made  in  this  research  are  discussed  in 
Chapter  III.  Chapter  IV  describes  the  results  from  the  simulations  of  seven  different 
scenarios.  Finally,  Chapter  V  wraps  everything  together  and  gives  ways  to  build  on 
the  results  of  this  thesis. 
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II.  Background 


2.1  Introduction 

The  goal  of  this  chapter  is  to  provide  the  reader  with  background  information 
about  orbits  in  space,  GPS,  and  some  existing  methods  for  either  using  GPS  to 
determine  position  or  reducing  errors  involved  with  GPS.  A  brief  overview  of  different 
orbits  will  be  described  with  GPS  reception  in  mind,  followed  by  an  overview  of  the 
current  GPS  constellation.  Afterwards,  errors  affecting  the  GPS  signal  are  reviewed 
before  presenting  the  position  estimation  method  and  error  sensitivity  method  to  be 
tested  with  this  research.  Finally,  other  known  methods  for  position  estimation  and 
error  reduction  will  be  given. 

2.2  Orbits 

Low  Earth  Orbit. 

Low  Earth  Orbit  (LEO)  is  a  stable  path  around  Earth  with  an  altitude  between 
80  km  and  2,000  km  [2],  The  orbital  velocity  for  satellites  in  LEO  average  around 
27,  500  km/hr  [17]  resulting  in  an  orbital  period  of  approximately  90  min.  LEO  con¬ 
tains  altitudes  that  are  within  the  GPS  antenna  main  beam’s  intended  area  of  effect, 
so  obtaining  multiple  GPS  satellites  in  view  for  position  determination  is  not  an  is¬ 
sue.  The  increased  velocity  compared  to  a  terrestrial  user  will  increase  the  range  of 
Doppler  values  seen  at  the  receiver.  Doppler  range  is  one  limiting  factor  in  choosing 
suitable  GPS  receivers. 

Higher  Orbits. 

There  are  two  orbits  at  altitudes  higher  than  LEO  that  are  of  interest  in  this 
research:  Medium  Earth  Orbit  (MEO)  and  GEO.  GEO  is  characterized  by  being  at 
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an  altitude  of  35,  786  km  and  having  an  orbital  period  equal  to  the  Earth’s  rotational 
period  [2],  The  MEO  falls  between  LEO  and  GEO,  and  has  an  altitude  between 
2,  000  km  and  35,  786  km  with  an  orbital  period  between  90  min  and  24hrs  [2], 

Satellites  operating  in  MEO  and  GEO  encounter  more  issues  with  GPS-based 
position  determination  than  those  in  LEO.  One  of  the  major  differences  in  comparison 
with  LEO  is  the  location  of  the  GPS  satellites  available  to  the  user.  For  a  user  located 
above  the  GPS  constellation,  the  GPS  signal  originates  from  beneath  the  user,  versus 
overhead  in  the  case  for  terrestrial  and  LEO  users  of  GPS  [18].  GPS  signals  that  are 
capable  of  reaching  the  user  will  originate  from  the  opposite  side  of  Earth  or  from 
the  antenna  side  lobes  of  nearby  GPS  SVs.  Figures  2.1  and  Figure  2.2  illustrate  the 
potential  GPS  transmit  beams  available  for  a  user  whose  orbital  altitude  is  above 
that  of  the  GPS  constellation. 
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Figure  2.2.  Potential  GPS  Beams  Available  at  HEO  [10] 


2.3  GPS 

Current  GPS  constellation  consists  of  31  operational  SVs  flying  in  MEO.  GPS 
SVs  fly  in  a  nearly  circular  orbit  with  an  altitude  of  approximately  20,  200  km  and 
an  orbital  period  of  approximately  12hrs  [14].  The  satellites  are  divided  into  six 
orbits,  each  containing  a  minimum  of  four  satellites.  Every  GPS  satellite  is  three-axis 
stabilized,  Earth  pointing,  with  an  inclination  of  55°.  The  system  is  now  mainly  used 
for  terrestrial  applications  but  has  recently  been  proven  to  work  with  receivers  located 
in  space.  GPS  receivers  are  being  flight  tested  on-board  certain  SVs  in  LEO,  e.g.  [7]. 

2.4  GPS  Error 

Signals  from  GPS  satellites  that  are  near  the  Earth’s  horizon,  though  unobstructed 
by  Earth,  still  have  errors  introduced  from  traversing  through  the  ionosphere.  The 
magnitude  of  errors  from  the  ionosphere  depends  on  the  distance  the  signal  travels 
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through  the  ionosphere  before  reaching  a  receiver.  A  good  illustration  of  this  is  from 
[18]  and  is  shown  in  Figure  2.3. 


GPS 


Figure  2.3.  Ionosphere  Layer  Illustration  (not  drawn  to  scale) 

Ionosphere  Effects. 

The  ionosphere  refers  to  a  region  of  ionized  gases  which  is  also  a  dispersive  medium 
for  radio  waves.  The  ionosphere  extends  from  around  50  km  to  1000  km  above 
the  Earth  [14].  When  the  GPS  signal  passes  through  the  ionosphere,  the  code  and 
modulating  signal  propagate  at  different  speeds.  This  results  in  the  code  phase  being 
delayed  and  the  carrier  phase  being  advanced  by  the  same  magnitude  [14] .  This  effect 
can  be  negated  by  using  a  dual  frequency  receiver  that  processes  both  GPS  LI  and 
L2  frequencies. 

For  the  purposes  of  determining  GPS  tracking  capability  above  LEO,  many  of 
the  current  research  makes  the  assumption  that  the  ionosphere  completely  obstructs 
GPS  signals  similar  to  the  Earth.  This  research  assumes  the  use  of  a  dual  frequency 
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receiver  so  ionospheric  effects  are  minimized.  More  information  on  the  method  can 
be  found  in  [14]. 

Signal  Degradation. 

As  a  signal  travels  through  a  space,  there  is  inherent  power  loss.  The  amount 
of  power  loss  is  directly  related  to  the  distance  traveled.  Received  signal  strength 
(RSS),  the  amount  of  signal  power  at  the  receiver,  from  the  it h  transmitter  can  be 
represented  as  follows  [21]: 

AM  =  Aq-  10np\ogw^~  +  ^  (2.1) 

r£o 

where  At(d,j)  is  the  estimate  of  received  power  in  Decibels  (dB)  at  a  distance  of  Rj,  R0 
is  the  reference  distance,  Ao  is  the  known  power  in  dB  at  the  reference  distance,  nr>  is 
the  path  loss  exponent,  and  Vt  is  a  zero- mean  Gaussian  random  variable  in  dB  [21]. 
GPS  uses  a  relatively  low  power  signal,  with  a  minimum  received  strength  for  a  3  dB 
polarized  antenna  on  Earth  between  — 160  dBW  and  — 153  dBW  for  the  LI  signal  [19]. 
In  [18],  one  of  the  underlying  assumptions  made  to  consider  an  unobstructed  GPS 
signal  available  to  the  simulated  GPS  receiver  at  GEO  is  for  the  power  received  to  be 
higher  than  —185  dBW.  GPS  signal  strength  is  referred  in  terms  of  the  Signal  to  Noise 
Ratio  (SNR)  or  more  commonly  Carrier  Power  to  Noise  Density  Ratio  (C/NO)  (in 
units  of  dBHz).  Comparing  the  signal  strength  to  the  background  noise  is  important 
when  determining  if  the  signal  is  strong  enough  to  be  acquired  by  a  receiver.  For 
a  simulation  study  in  [13],  a  usable  GPS  signal  required  the  C/NO  at  the  receiver 
to  remain  above  the  user  defined  minimum  long  enough  to  decode  the  navigation 
message.  Another  example  is  an  actual  flight  test  in  [4]  where  the  GPS  receiver  is 
equipped  on  a  satellite  named  Equator-S.  The  receiver  had  an  acquisition  threshold 
of  37  dBHz  but  did  not  directly  measure  and  output  SNR,  so  a  relation  between  its 


signal  strength  parameter  and  SNR  value  was  used  to  determine  relative  physical  SNR 
values.  This  helps  portray  the  importance  of  keeping  track  of  SNR  and  selecting  a 
receiver  with  an  appropriate  threshold. 

2.5  Position  Determination 

Navigation  data  in  the  GPS  signal  contains  a  variety  of  information  including 
signal  transmission  times,  clock  error  corrections,  ephemeris  data  for  that  specific 
GPS  SV,  almanac  information  for  the  entire  GPS  constellation,  satellite  state  of 
health,  etc.  Signal  transmission  time,  clock  error  corrections,  and  ephemeris  data  are 
used  to  calculate  the  location  of  each  GPS  SV  in  view  as  well  as  approximate  ranges 
from  the  SVs  to  the  user. 

Approximate  range  measurements  derived  from  the  difference  between  GPS  signal 
transmit  time  and  time  of  arrival  at  the  receiver  are  called  pseudoranges.  A  pseu¬ 
dorange  is  an  approximate  range  from  the  user’s  current  position  to  the  GPS  SV’s 
position.  By  using  the  pseudoranges  from  multiple  SVs  together  with  knowledge  of 
SV  location,  a  receiver  can  estimate  its  position. 

Multiple  algorithms  are  openly  available  that  use  pseudoranges  to  estimate  a  user’s 
position.  One  position  determination  algorithm  in  particular,  generally  referred  to 
as  the  Newton- Raphson  method,  requires  a  minimum  of  four  GPS  SVs  in  view.  The 
Newton-Raphson  method  is  implemented  in  this  research  and  a  summary  of  this 
method  is  given  in  the  next  section  based  on  a  more  in-depth  description  in  [14]. 

Newton  Raphson  -  Least  Squares  Algorithm. 

The  Newton-Raphson  method  is  based  around  solving  the  range  equation.  The 
process  of  estimating  user  position  involves  calculating  range  to  each  GPS  SV  based 
on  a  nominal  initial  guess  of  the  current  user’s  position.  The  calculated  range  is 
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compared  with  the  measured  pseudorange  to  produce  an  error.  Afterwards,  the  first 
order  Taylor  series  expansion  of  the  range  equation  in  Equation  (2.2)  is  used  to  show 
sensitivities  to  a  user’s  position.  The  last  step  combines  the  sensitivity  information 
with  at  least  four  of  the  calculated  errors  to  solve  for  a  minimum  least  squares  position 
error. 

The  Newton- Raphson  method  requires  a  minimum  of  four  satellites  because  it 
involves  solving  a  linearized  form  of  the  ranging  equation: 

P[u]  =  \! (TW  -  £u)2  +  (p[s]  -  Vu)2  +  (^[s]  -  ~Zu)2  +  C*tu  (2.2) 

about  an  approximate  user’s  position  and  receiver  clock  error.  In  Equation  (2.2)  the 
superscript  (s)  represents  the  GPS  SV,  the  subscript  u  refers  to  the  user,  p  is  the 
range,  x,y,  z  are  position  axes  and  c*tu  is  the  receiver  clock  error  tu  being  multiplied 
by  the  speed  of  light  c  to  convert  it  into  units  of  length.  After  the  range  equation  is 
linearized,  it  is  iteratively  solved  for  the  four  unknowns  (user’s  three  dimensions  of 
position  and  receiver  clock  error)  [14], 

An  initial  estimate  of  user  position  [xu,yu,zu\  and  receiver  clock  error  tu  are  re¬ 
quired  for  the  Newton-Raphson  method  and  will  be  denoted  in  equations  with  a 
subscript  of  0.  In  the  updated  range  equation: 

Po 1  =  \J (^[s]  -  £o)2  +  (p[s]  -  Po)2  +  (z[s]  -  A))2  +  C*t0  (2.3) 

Pq- '  represents  the  range  from  the  sth  GPS  SV  to  the  initial  estimated  user  position  at 
[xo,  Ho,  zo]  affected  by  an  estimated  receiver  clock  error  of  to-  Relating  the  estimated 
values  to  the  truth  is  accomplished  by  using  the  variables  Sx  and  St  to  represent  the 
difference  between  truth  and  estimation.  Expressing  the  true  position  on  each  axis  as 
xu  =  x0  +  Sx,  yu  =  y0  +  Sy,  zu  =  z0  +  Sz  and  true  receiver  clock  error  as  /„  =  t0  f  St 
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allows  Equation  (2.3)  to  be  rewritten  in  terms  the  error  St  and  Sx  as  shown  below: 


Pu]  =  V  (x[s]  -  x0  -  Sx)2  +  (yb]  —  |/o  —  Sy)2  +  (z^  -  z0  -  Sz)2  +  c(tQ  +  St)  (2.4) 


The  next  step  is  to  linearize  Equation  (2.4)  by  taking  the  Taylor  series  expansion 
of  it  and  only  using  the  first  order  terms.  The  resulting  form  is  shown  in  Equation 
(2.5)  below: 


pW  =  P?  -  A*  -  y—TT^  Ap  -  +  cAt 


where  Ax,  Ay,  A z,  and  At  represent  the  approximated  difference  between  the  initial 


estimated  user  position  and  the  truth.  Expressing  the  range  of  all  GPS  SVs  in  view 
in  the  form  of  Equation  (2.5)  allows  the  relationship  to  be  expressed  in  the  form  of  a 
matrix  operation  with  some  rearranging.  The  resulting  matrix  operation  is  shown  in 
Equation  (2.6)  and  an  expanded  form  is  shown  in  Equation  (2.7)  [20]. 


A  p  =  HAx  (2.6) 

Apt1!  —1  Ax 

Apt2!  ajf  a $  —1  Ay 

\  \  A  z 

Ap^J  |_af]  of1  af1  -lj  [cAt 

Ap  is  a  vector  containing  the  differences  between  estimated  ranges  and  the  true 
ranges  pf^  —  p^ ,  H  is  a  matrix  made  up  of  the  Taylor  series  expansion  for  each  GPS 
SV,  Ax  is  a  vector  representing  the  difference  between  the  true  user  position  and 
the  estimated  position,  and  N  is  the  number  of  visible  satellites.  The  Taylor  series 
components  represented  in  the  H  matrix  are  shown  in  Equation  (2.8). 
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aH  = 


aW  -  y- 

ay 


Vo 


aW  = 


(2-8) 

Po  J  Po  J  Po  J 

Since  Ax  is  the  difference  between  the  true  user  position  and  the  estimated  po¬ 
sition,  it  can  be  solved  for  and  used  to  correct  the  initial  position  estimate  using 
Equations  (2.9)  and  (2.10): 


Ax  =  (HrH)_1HrAp  (2.9) 

x  =  xo  +  Ax  (2-10) 

where  x  is  the  new  estimate.  Simple  matrix  operations  such  as  this  can  be  applied  due 
to  the  use  of  linearized  estimates  and  neglected  higher  order  terms.  The  whole  process 
is  now  required  to  be  iterated  for  those  same  reasons.  The  iterative  process  involves 
using  an  initial  position  estimate,  solving  for  Ax,  updating  the  position  estimate  x, 
using  the  new  estimate  to  start  the  process  over  again,  and  repeating  until  the  least 
squares  of  Ax  is  under  a  set  threshold. 

2.6  Geometry  Errors 

Precision  of  the  Least  Squares  method  is  not  only  dependent  on  the  number  of 
available  GPS  SVs  but  also  the  geometrical  orientation  of  those  SVs  relative  to  the 
user.  A  position  and  time  solution  cannot  be  calculated  if  less  than  four  GPS  SVs 
are  available.  Even  with  four  or  more  GPS  SVs  present,  position  precision  can  be 
degraded  due  to  poor  geometry  if  all  the  SVs  are  in  the  same  relative  direction.  A 
two  dimensional  example  of  this  is  presented  in  Figure  2.4.  DOP  is  a  metric  that 
takes  into  account  the  effect  of  user-satellite  geometry  on  errors.  The  more  favorable 
the  user-satellite  geometry  is,  the  lower  the  DOP  value  [14]. 
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Good  Geometry 


Bad  Geometry 


S  =  User 
"  =  Range  Error 
“  =  Possible  Solutions 


Figure  2.4.  Illustration  of  the  effect  of  geometry  on  position  error 

Dilution  of  Precision. 

DOP  is  a  unitless  measurement  that  directly  relates  sensitivity  of  range  measure¬ 
ment  errors  to  user  position  errors.  Range  measurements  in  this  context  refer  to  the 
pseudoranges  from  user  to  GPS  SVs.  The  relation  between  range  and  position  errors 
takes  into  account  the  geometry  of  the  GPS  SVs  relative  to  the  user’s  position.  DOP 
measurement  favors  relatively  large  angles  of  separation  between  GPS  SVs  visible  to 
the  user.  DOP  originates  as  a  4x4  covariance  matrix  relating  to  the  four  dimensions 
of  position  (including  time) .  This  matrix  is  calculated  in  the  process  of  solving  for  Ax 
in  the  Newton  Raphson  iterative  process.  In  Equation  (2.9),  the  (HTH)-1  calculation 
produces  the  DOP  matrix  as  shown  in  Equation  (2.11). 
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Dxx 

DXy 

Dxz 

Dxt 

DOP  Matrix 

=  (HtH)-1  = 

Dyx 

Dyy 

Dyz 

Dyt 

(2.11) 

DZX 

DZy 

Dzz 

Dzt 

Dtx 

Dty 

Dtz 

Dtt  _ 

Each  diagonal  represents  the 

variance  of  one 

axis, 

Dxx  — 

z 

Dyy  = 

al 

etc.  It  is  not 

practical  to  refer  to  the  DOP  matrix  as  a  whole  so  DOP  is  normally  referred  to  with  a 
scalar  value.  The  scalar  value  is  derived  from  a  combination  of  diagonal  values  in  the 
4x4  matrix.  Two  common  combinations  are  Geometric  Dilution  of  Precision  (GDOP), 
which  includes  all  four  dimensions;  and  Position  Dilution  of  Precision  (PDOP),  which 
includes  all  but  time.  The  equations  for  the  different  types  of  DOP  are  [20] : 

Geometric  DOP  =  yf . Dxx  +  ~D^y  +  Dzz  +  Dtt 

Position  DOP  =  yf Dxx  +  ~D^y  +  Dzz 

Horizontal  DOP  =  yj Dxx  +  _DW  (2.12) 

Vertical  DOP  =  y/ Dzz 

Time  DOP  =  yj Dtt 

As  a  note,  the  Horizontal  and  Vertical  DOP  is  only  applicable  when  using  the  East- 
North-Up  (ENU)  or  a  similar  coordinate  system.  The  lower  the  DOP  value  the  better 
the  positional  precision,  with  a  DOP  of  1.0  being  ideal.  DOP  should  be  viewed  as  a 
multiplier  for  the  potential  user  position  error.  For  example  if  the  user  had  a  range 
error  of  5  m  then  a  GDOP  of  5  would  indicate  the  potential  for  up  to  25  m  error  in 
the  user’s  position. 
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2.7  Previous  Research  Incorporating  an  On-Board  GPS  Receiver 

A  GPS  receiver  was  installed  on-board  the  Highly  Elliptical  Orbit  (HEO)  satel¬ 
lite  Equator-S  and  launched  in  December  1998  into  a  Geostationary  Transfer  Or¬ 
bit  (GTO)  [4],  The  main  goal  of  the  experiment  was  to  determine  if  reception  of 
GPS  signal  was  possible  above  the  GPS  constellation.  Despite  various  configuration 
issues  that  arose,  the  research  proved  that  tracking  GPS  at  altitudes  above  the  GPS 
constellation  and  near  GEO  is  possible  [4],  The  research  also  proved  that  reception 
of  signals  originating  from  the  side  lobes  of  the  GPS  transmit  antenna  is  also  possible 

[4]- 

The  US  Air  Force  Academy  sponsored  an  experiment  called  “Falcon  Gold”  to 
measure  GPS  signals  at  high  orbital  altitudes  using  the  NAVSYS  TIDGET  sensor 
[16].  Falcon  Gold  was  flown  in  November  1997,  mounted  to  the  side  of  the  Centaur 
upper  stage.  Once  the  Centaur  upper  stage  separated  from  the  payload  in  GTO,  the 
Falcon  Gold  sensor  started  collecting  40  ms  of  sampled  data  every  5  minutes.  The 
data  was  not  processed  onboard.  This  experiment  confirmed  that  a  low  cost  sensor 
can  be  used  to  detect  GPS  signals  in  GTO,  including  those  signals  from  the  GPS 
antenna  side  lobes.  The  research  also  provides  the  C/NO  values  observed  at  various 
altitudes  in  the  GTO  orbit. 

There  are  other  research  that  use  a  GPS  receiver  on-board  a  satellite  as  part  of 
a  bigger  experiment.  The  “Navigation  and  Occultation  experiment”  uses  a  GPS 
receiver  on  a  small  satellite  “TET-1”  to  prove  that  low  cost  commercial-off-the-shelf 
hardware  can  be  used  in  space-borne  applications  with  only  minor  changes  to  the 
receiver’s  firmware  [7].  While  the  GPS  receivers  on  the  “Constellation  Observing 
System  for  Meteorology,  Ionosphere  and  Climate  (COSMIC)”  satellites  are  used  for 
precise  orbit  determination  and  surveying  the  environment  through  radio  occultation 
[12]- 
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2.8  Tools  Used  in  Research 


A  combination  of  three  different  tools  are  used  for  this  research.  These  tools 
include  Systems  Tool  Kit  (STK),  used  for  orbit  propagation;  MATLAB,  for  numerical 
computation  and  object  oriented  programming;  and  the  SPIRENT  GSS8000,  for 
Radio  Frequency  (RF)  signal  generation.  Each  of  these  systems  are  described  in 
more  detail  below. 

STK. 

STK  is  a  modeling  and  simulation  program  from  Analytical  Graphics  Incorpo¬ 
rated  (AGI),  located  in  Exton  PA,  that  deals  with  aircraft,  satellite,  ground  vehicles, 
sensors,  as  well  as  communication  elements.  STK  is  known  for  its  ability  to  model 
orbit  propagations.  The  program  offers  a  variety  of  different  orbit  propagators  that 
differ  in  fidelity  and  complexity  of  included  orbital  perturbations.  Using  STK  allows 
for  study  of  orbit  models  and  parameter  sensitivity.  Loading  a  satellite,  the  GPS 
constellation,  or  an  antenna  into  STK  can  be  done  from  an  external  file  or  STK’s  de¬ 
fault  repository.  There  are  also  options  for  calculating  connection  times,  dilution  of 
precision,  navigation  accuracy,  as  well  as  setting  various  constraints.  Since  these  func¬ 
tions  are  used  in  industry  for  terrestrial  scenarios  involving  GPS,  STK  will  be  used 
to  validate  the  software  developed  for  this  research  using  similar  terrestrial  scenarios. 

MATLAB. 

MATLAB  is  a  high  level  language  and  numerical  computing  environment.  MAT¬ 
LAB  does  not  have  a  built  in  orbit  propagator  but  it  does  have  the  ability  to  integrate 
with  STK.  STK  and  MATLAB  can  be  integrated  through  file  interoperability,  plug¬ 
ins,  and  through  Microsoft’s  COM  library  which  allows  MATLAB  to  operate  STK’s 
functions  using  object-oriented  programming  commands.  For  this  research,  MAT- 
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LAB  is  used  to  automate  the  scenario  creation  in  STK,  load  and  propagate  the  user 
and  GPS  orbits  at  the  desired  time  step,  attach  transmit  and  receive  sensors  to  each 
satellite,  and  finally  export  data  back  into  MATLAB  for  processing.  This  level  of  au¬ 
tomation  to  do  the  systematic  sweeps  needed  for  parameter  studies  allows  this  study 
to  be  feasible  based  on  available  resources. 


MATLAB  GUI. 


Using  MATLAB  to  connect  to  STK  with  connect  commands  involves  a  moderate 
learning  curve,  partly  due  to  the  complexity  in  deciphering  errors  and  the  translation 
of  examples  given  in  C#  into  MATLAB’s  COM  syntax.  AFIT  Summer  2014  interns 
Arnett,  Wagner,  and  Dieter le  created  a  Graphic  User  Interface  (GUI)  to  handle  the 
MATLAB  to  STK  transactions.  Figure  2.5  shows  the  top  level  interface  for  the  GUI. 


I  a  I  alB 


4gHT 


Target  Satellite  Inputs 

Semi-Major  Axis  (km)  |  7096 
Eccentricity  p 

Inclination  (deg)  24 

Argument  of  Perigee  (deg)  j  p 
RAAN(deg)  p 

True  Anomaly  (deg)  p 


Additional  Options 


GNSS  Constellations  TLE/Almanac  Files 

[7]  GPS 

C:\Users\Anthony\DesktopVGUkyuma16May.aim  [77] 

0  GLONASS 

C:\Users\Anthony\DesktopVGUtglonass.tce  |  [77] 

0  Beidou 

C:\Users\Anthony\Desktop\GUtoeidou.tce  |  [77] 

0  Galileo 

C:\Users\Anthony\Desktop\GUHgalileo.tce  |  [77] 

0  Calculate  and  store  GNSS  positions 

Simulation  Time 


Generate  Reports 

3  Sim  □  UMT  0  GGA  0  ZDA 


Generate  Plots 

0 


Simulation  Time  umestep  [  60  (sec) 

Start  Date  P9Jun2014  (<M  mmm  yyyy)  Start  Time  12:00:00  (hh:mm:ss) 

End  Date  09Jun2014  (<M  mmm  yyyy)  End  Time  12:24:00  (hh:mm:ss) 


Plate  Carree  Ground  Track 


Figure  2.5.  MATLAB  to  STK  GUI 


The  GUI  requires  the  user  to  input  the  satellite’s  six  classical  orbital  elements, 
simulation  time  and  time  step,  GPS  almanac  location  or  other  Global  Navigation 
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Satellite  System  (GNSS)  Two-Line  Element  (TLE)  file  locations,  and  to  select  which 
reports  are  to  be  generated.  Some  of  the  default  parameters  that  can  be  changed 
in  the  GUI  include  the  receiver  antenna’s  pattern,  receiver  antenna  orientation  and 
offset,  C/NO  constraint,  and  orbit  propagator  preference. 

SPIRENT  GSS8000. 

SPIRENT  GSS8000  is  a  RF  constellation  simulator  that  supports  GPS,  Galileo, 
Global’naya  Navigatsionnaya  Sputnikovaya  Sistema  (GLONASS),  and  Quasi-Zenith 
Satellite  System  GNSSs.  It  is  designed  as  a  platform  for  the  development  and  testing 
of  GNSS  enabled  devices  as  well  as  the  GNSS  technology  itself.  Spirent  Communica¬ 
tions  has  developed  different  models  of  GNSS  signal  simulators,  each  with  different 
levels  of  functionality.  Utilization  of  SPIRENT’s  RF  generator,  as  well  as  various 
models  included  in  the  SPIRENT  software  package,  is  found  in  different  research  ef¬ 
forts  such  as  in  [5]  and  [6].  The  simulator  is  able  to  produce  the  GNSS  RF  signal 
expected  to  be  received  at  the  user’s  location  by  taking  into  account  variables  such 
as:  user  position  and  velocity,  antenna  patterns,  obscuration  effects,  signal  propaga¬ 
tion,  and  various  error  models  (including  multipath).  The  RF  signal  is  generated  in 
real  time  and  transmitted  via  antenna  or  cable  to  a  GNSS  receiver  for  processing. 
A  GNSS  enabled  device  can  process  the  signal  in  real  time  as  the  signal  is  being 
transmitted.  The  signal  can  also  be  recorded  and  downsampled  with  an  RF  front  end 
for  eventual  post-processing  using  a  Software  Defined  Radio  (SDR).  The  simulator’s 
real  time  signal  generation  requires  the  time  step  for  the  user  motion  input  to  be  in 
one  second  increments  or  less.  If  using  one  second  increments,  an  option  can  be  set 
to  allow  the  simulator  to  interpolate  the  user’s  position  in  the  time  spans  between 
the  one  second  increments. 
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GPS  Antenna  Model. 


SPIRENT  specializes  in  RF  signal  simulator  so  the  simulator  includes  a  higher 
fidelity  model  of  the  GPS  transmit  antenna  beams.  Unlike  STK,  which  only  models 
the  main  beam  of  the  GPS  antenna,  SPIRENT’s  model  includes  the  antenna  sidelobes 
as  well.  STK’s  antenna  model  is  included  by  default  with  their  GPS  transmitter 
model,  as  shown  in  Figure  2.6. 


Figure  2.6.  STK’s  GPS  Transmit  Antenna  Beam  Model 

SPIRENT’s  antenna  model  offers  a  more  detailed  version  of  the  GPS  transmit 
beam  with  a  variety  of  attenuation  values.  Figure  2.7  shows  a  visualization  of  the 
GPS  transmit  antenna’s  attenuation  model  used  in  SPIRENT.  The  model  is  also 
exportable  as  a  table  of  attenuation  values  by  azimuth  and  elevation  which  is  more 
useful  for  increasing  simulation  fidelity.  For  the  purposes  of  this  research,  SPIRENTs 
antenna  model  is  imported  into  MATLAB  and  used  to  calculate  the  GPS  signal 
power  expected  at  the  receiver.  The  antenna  pattern  output  from  SPIRENT’s  model 
resembles  the  GPS  signal  model  shown  in  Figure  2.2. 
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Figure  2.7.  SPIRENT’s  GPS  Transmit  Antenna  Beam  Model 


2.9  Summary 

The  chapter  finished  with  the  tools  that  will  be  used  to  simulate  the  GPS  -  user 
interactions  that  will  result  in  evaluating  the  position  estimation  errors.  Various  po¬ 
sition  estimators  are  currently  available  in  the  PNT  community  and  a  few  common 
techniques  were  mentioned.  This  research  chose  to  implement  the  Newton-Raphson 
position  estimation  algorithm  in  order  to  evaluate  both  position  errors  and  the  appli¬ 
cability  of  the  DOP  metric.  DOP  is  a  metric  that  evaluates  position  precision  errors 
based  on  the  geometry  of  GPS  SVs  relative  to  the  user.  Space  based  errors  affecting 
the  GPS  signal  were  also  discussed  along  with  a  brief  description  of  the  orbits  in  space 
of  interest  to  this  research. 
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III.  Methodology 


3.1  Introduction 

Before  position  can  be  calculated,  it  is  necessary  to  determine  the  number  of 
GPS  SVs  the  receiver  is  able  to  detect.  To  be  detectable,  the  signal  from  the  GPS 
to  the  receiver  must  not  be  obstructed  by  Earth,  fall  within  the  receiver’s  Doppler 
band  limits,  and  be  at  a  sufficient  power  level  at  time  of  reception.  This  chapter  will 
discuss  the  initial  setup  in  STK,  line  of  sight  vector  calculation,  Doppler  and  power 
calculation,  and  finally  position  calculations. 

3.2  Work  Flow 

Figure  3.1  shows  the  flow  of  this  research  from  start  to  finish.  The  scenario  specific 
data  is  inputted  into  the  MATLAB  to  STK  GUI.  This  GUI  is  responsible  for  setting 
up  and  running  the  scenario  in  STK.  STK  is  a  closed  source  commercial  product 
so  the  source  code  is  not  released  to  the  public.  As  a  result,  this  research  only  uses 
STK  to  1)  propagate  the  orbits  of  the  user  and  GPS  satellites,  and  2)  return  the 
Earth  Centered  Earth  Fixed  (ECEF)  or  Earth  Centered  Inertial  (ECI)  position  and 
velocities  for  each  time  step.  Next,  the  GUI  takes  the  outputs  from  STK,  formats 
them  in  MATLAB,  and  then  saves  off  two  separate  data  files.  The  position  and 
velocity  file  is  used  for  analysis  in  MATLAB  as  described  in  Section  3.5,  while  the 
user  motion  trajectory  file  is  used  for  physical  simulation  of  the  RF  signal  with  the 
GNSS  simulator  in  Section  3.8.  Both  options  result  in  pseudoranges  to  be  used  in  the 
Newton- Raphson  least  squares  position  calculation  algorithm.  To  determine  position 
error,  the  results  of  the  position  calculation  algorithm  is  compared  with  the  true  user 
position  from  STK. 
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Figure  3.1.  Research  Work  Flow 


3.3  Assumptions 

A  few  assumptions  are  made  when  building  the  scenarios  to  narrow  the  scope  of 
the  research  or  to  ensure  realistic  values  are  achieved  for  noise  and  equipment  settings. 
Below  is  a  list  of  the  assumptions. 

1.  The  receiver  antenna  is  isotropic,  thus  able  to  receive  signals  from  360°. 

2.  Range  measurements  are  not  affected  by  multipath  errors. 

3.  Dual  frequency  (LI  and  L2  bands)  measurements  are  used  in  the  receiver  to 
implement  ionosphere  free  measurements.  Outside  of  ionosphere  free  measure¬ 
ments,  the  receiver  will  only  process  the  LI  civilian  frequency. 
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4.  The  Earth’s  radius  will  include  the  height  of  the  troposphere. 


5.  A  third  order  Phase  Lock  Loop  (PLL)  is  used  to  accommodate  the  high  dy¬ 
namics  of  a  elliptical  orbit. 

3.4  Code  Validation 

Coding  for  this  research  takes  place  in  MATLAB  in  order  to  provide  a  clear  view 
of  all  calculations  taking  place.  All  calculations  are  manually  coded  from  known 
equations  in  GNSS  open  literature,  such  as  Line  of  Sight  (LOS),  DOP,  and  Doppler. 
Confidence  in  the  MATLAB  code  working  as  expected  was  gained  by  comparing  re¬ 
sults  with  STK  for  a  scenario  involving  terrestrial  based  receiver  and  no  consideration 
for  GPS  antenna  side  lobes.  Even  though  STK  is  only  being  used  for  the  position 
and  velocity  information,  it  is  capable  of  providing  information  on  DOP,  C/NO,  nav¬ 
igation  accuracy,  and  transmitter  to  receiver  connection  time.  Different  analyses  on 
the  data  can  be  run  in  STK  by  using  the  Reports  &  Graph  manager.  These  values 
can  then  be  compared  with  the  results  from  MATLAB  to  verify  similarity  of  both 
processes. 

The  following  sections  provide  more  detail  about  what  happens  during  the  sup¬ 
plementary  analysis  phase  and  the  physical  simulation  phase.  A  detailed  flow  chart 
of  the  supplementary  analysis  is  shown  in  Figure  3.2. 
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Figure  3.2.  Detailed  Supplementary  Analysis 


3.5  Supplementary  Analysis 

Supplementary  analysis  occurs  in  MATLAB  and  handles  all  of  the  calculations, 
criteria  comparison,  position  estimation,  and  finally  data  analysis  involved  in  this 
research.  Required  inputs  are  position  and  velocity  vectors  from  STK,  GPS  antenna 
model,  and  the  GPS  almanac  file  associated  with  the  scenario  time  frame.  The  first 
portion  of  the  supplementary  analysis  is  to  determine  the  GPS  satellites  to  which 
the  user  connects.  A  connection  is  considered  established  if  the  user  has  LOS  to  the 
GPS  SV,  a  Doppler  shift  magnitude  under  a  predefined  threshold,  and  received  signal 
strength  above  a  predefined  level. 
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Line  of  Sight. 


For  two  objects  to  have  a  Line  of  Sight  connection  in  this  research  the  vector 
connecting  both  objects  is  required  to  be  unobstructed  by  Earth.  The  connecting  LOS 
vector  Plos  is  obtained  by  subtracting  the  user’s  position  vector  Py  from  the  selected 
GPS  SV’s  vector  Pgps-  Figure  3.3  illustrates  this  operation:  Plos  =  Pu  —  Pgps- 


>  User  Position 
GPS  Position 
LOS  Vector 

Figure  3.3.  LOS  Vector  Calculation 


Next,  the  range  of  the  shortest  distance  vector  from  Earth  to  the  LOS  vector  is  ex¬ 
amined.  The  shortest  distance  vector  is  found  by  first  calculating  the  projection 
of  the  user  position  vector  onto  the  LOS  unit  vector  Plos  using  the  dot  prod¬ 
uct.  Afterwards,  that  projection  vector  is  subtracted  from  the  original  user  posi¬ 
tion  vector  to  obtain  the  shortest  distance  vector.  The  range  is  found  by  taking  the 
norm  of  the  shortest  distance  vector.  An  illustration  of  this  mathematical  operation, 
REarth2Signal  =  ||P[/  -  ((Plos  •  Pp)Plos)||,  is  shown  in  Figure  3.4. 
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User  projection  on  the  LOS  vector 


User  Position 
GPS  Position 
LOS  Vector 
Projection  Vector 
Earth  to  LOS  Vector 


User  position  vector  -  Projection  vector 
=  Earth  to  LOS  vector 


Figure  3.4.  Earth  to  LOS  Vector  Calculation 


Normally  the  projection  vector  is  calculated  with  a  dot  product.  However,  for  this 
research  the  dot  product  to  find  the  user-to-LOS  projection  is  limited  to  account  for 
two  outlier  cases,  when  the  dot  product  is  negative  or  the  magnitude  is  greater  than 
the  LOS  magnitude. 

A  negative  dot  product  means  the  projection  of  the  LOS  vector  onto  the  user 
position  vector  is  pointing  towards  the  Earth.  This  happens  when  the  LOS  vector’s 
origin,  i.e.  GPS  SV’s  position,  is  at  a  higher  altitude  relative  to  Earth  than  the 
user.  If  the  GPS  SV  is  above  the  user,  relative  to  Earth,  then  obstruction  from  Earth 
cannot  occur.  Thus  the  first  limitation  is  to  set  a  lower  bound  of  zero  on  the  dot 
product  by  taking  the  maximum  value  between  zero  and  the  calculated  dot  product 
to  ensure  non-negative  values.  In  other  words,  the  zero  lower  bound  prevents  the 
calculated  nearest  point  on  the  LOS  vector  from  being  beyond  the  user  satellite. 

The  second  limitation  involves  setting  an  upper  bound  for  the  projection  vector. 
This  upper  bound  prevents  the  nearest  point  on  the  LOS  vector  from  being  beyond 
the  GPS  SV.  If  the  magnitude  of  the  LOS  vector  is  smaller  than  the  magnitude 
of  the  user’s  position  vector,  then  the  possibility  exists  for  the  projection  vector  to 
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be  greater  than  the  LOS  vector.  This  in  turn  causes  problems  when  calculating  the 
shortest  distance  vector  from  Earth  to  the  LOS  vector.  Figure  3.5  illustrates  a  simple 
example  of  the  reason  for  the  projection  vector  upper  bound.  In  this  example,  Earth, 
the  GPS  SV,  and  the  user  are  all  inline  with  the  user  being  above  the  GPS  SV. 


Earth  to  LOS  vector  is 
same  magnitude  as 
GPS  position 


User  Position  Vector 
GPS  Position  Vector 
LOS  Vector 
Projection  Vector 
Earth  to  LOS  Vector 


Figure  3.5.  Unbounded  vs  Bounded  projection  comparison 


In  this  situation  the  magnitude  of  the  user  position  vector  is  greater  than  that  of  the 
LOS  vector.  When  the  user  position-to-LOS  projection  is  found  and  subtracted  from 
the  user  position  vector,  the  resulting  Earth  to  LOS  vector  will  point  to  a  location 
outside  of  the  bounds  of  the  LOS  end  points.  Once  the  Earth-to-LOS  vector  is  found, 
the  magnitude  is  used  to  determine  if  a  GPS  SV  is  considered  in  view  to  the  user. 
This  will  be  discussed  in  more  detail  in  Section  3.6. 
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Doppler  Shift. 


The  receiver’s  incoming  GPS  RF  signal  is  a  combined  signal  comprised  of  trans¬ 
missions  from  multiple  GPS  SVs.  To  determine  which  GPS  signals  are  present,  the 
receiver  produces  a  reference  signal  for  each  GPS  SV  to  correlate  with  the  incoming 
RF  signal.  The  receiver  has  to  account  for  the  change  in  frequency  due  to  Doppler 
effects  from  the  relative  speed,  so  the  reference  signal’s  frequency  is  varied  over  a 
range  of  potential  Doppler  values.  If  the  incoming  signal  has  a  Doppler  shift  outside 
of  range  searched,  then  the  signal  is  not  recognized  by  the  receiver.  Thus  for  this 
research,  the  Doppler  effect  caused  by  relative  velocity  is  calculated  for  each  GPS 
SV. 

Doppler  shift  depends  on  the  velocity  at  which  two  objects  are  approaching  or 
departing  one  another.  This  value  is  the  same  as  the  velocity  along  the  LOS  vector. 
First,  the  relative  velocity  for  both  the  user  and  GPS  SV  is  calculated  by  subtracting 
the  velocity  vectors.  The  relative  velocity  vector  is  then  projected  onto  the  LOS 
vector  to  determine  the  magnitude  of  approach  or  separation. 


GPS  Velocity  vector  -  User  Velocity  vector 
=  Relative  Velocity 


User  Velocity  Vector 
GPS  Velocity  Vector 
LOS  Vector 
Relative  Velocity 

Projected  Relative  Velocity  Vector 


Relative  Velocity  vector  projected  onto 
the  LOS  vector 


Figure  3.6.  Relative  Velocity  Calculation 


Doppler  shift  is  calculated  as  follows: 
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Doppler 


(3.1) 


VlosII  *  Jli 
c 

where  V los  is  the  projected  velocity  on  the  LOS  vector,  fn  is  the  frequency  of  the 
LI  band:  1575.42  MHz,  and  c  is  the  speed  of  light.  The  Doppler  shift  value  is  then 
compared  against  a  set  criteria  to  determine  if  the  receiver  is  able  to  process  the 
signal  from  the  GPS  SV.  More  details  on  this  process  can  be  found  in  Section  3.6. 

Received  Signal  Strength. 

Another  case  where  a  GPS  signal  is  not  seen  by  the  receiver  is  due  to  low  signal 
strength.  If  the  received  GPS  signal  strength  is  too  low  then  it  will  be  indistinguish¬ 
able  from  the  background  noise  even  after  signal  processing.  Adequate  signal  strength 
in  this  research  is  determined  in  two  different  ways,  depending  on  whether  the  GPS 
antenna  sidelobes  are  being  considered,  or  only  the  main  beam. 

No  Sidelobes. 

The  simplest  case  considered  is  when  only  the  main  beam  is  used  (so  no  sidelobes). 
In  these  scenarios,  the  received  GPS  signal  is  considered  strong  enough  to  be  processed 
if  the  user  is  located  within  the  main  beam  path  of  that  GPS  SV.  Equation  (3.2) 
shows  the  calculation  of  the  GPS  to  user  angle. 

@G2U  =  COS  X(  —  PgPS  •  Flos)  (2-2) 

The  angle  represents  how  far  the  user  is  from  the  GPS  antenna’s  boresight,  which  can 
be  compared  to  the  main  beam  half  cone  angle  listed  in  the  GPS  Interface  Control 
Document  (ICD)  [19]. 
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Sidelobes  Considered. 


Power  in  the  sidelobes  is  reduced  when  compared  to  the  output  from  the  main 
beam  but  can  still  be  strong  enough  for  signal  reception.  When  considering  the 
effects  of  sidelobes,  a  different  approach  for  determining  adequate  signal  strength  is 
required.  Instead  of  having  the  boresight  angle  determine  reception,  the  boresight 
angle  determines  the  amount  of  gain/attenuation  applied  from  the  GPS  antenna  in 
dB.  The  gain  values  are  generated  by  interpolating  the  GPS  antenna  attenuation  table 
exported  from  SPIRENT.  Normally  the  next  step  would  be  to  input  the  antenna  gain 
and  signal  transmission  power  into  a  free-space  path  loss  equation  to  determine  signal 
strength  at  reception.  Unfortunately,  the  GPS  ICD  does  not  mention  specifics  on  the 
GPS  signal  transmission  power,  only  the  lowest  guaranteed  power  at  reception  on 
Earth  [19].  Thus  a  modified  path  loss  equation  from  SPIRENT  [1]  is  implemented  to 
determine  the  signal  power  at  reception,  shown  in  Equation  (3.3). 

SPj-cot-  =  SP ICD  +  Og  +  20  *  log10(  — - )  +  Gtx  +  Grx  (3-3) 

J^los 

In  this  equation,  SPrcm.  is  the  power  at  the  receiver  in  dBm,  SP icd  is  the  lowest 
guaranteed  GPS  power  at  reception  mentioned  in  the  GPS  ICD,  Og  is  a  global  offset 
constant  used  in  SPIRENT,  R0  is  the  range  of  a  GPS  SV  at  the  furthest  point  away 
from  a  terrestrial  user,  Rlos  is  the  actual  range  from  the  GPS  SV  to  the  user,  Gtx 
is  the  gain  from  the  transmission  antenna,  and  Grx  is  the  receiver  antenna  gain  [1]. 

SP icd  value  is  fisted  as  —128.5  dBm  [19],  the  Og  default  value  in  SPIRENT  is 
10  dB  [1],  Rlos  is  the  magnitude  of  Flos',  Gtx  is  the  value  from  comparing  GPS 
boresight  angle  with  the  antenna  model  interpolation,  and  the  receiver  is  assumed  to 
be  isotropic  resulting  in  Grx  being  set  to  0  dB.  When  calculating  R0  it  is  important 
to  note  that  the  farthest  visible  GPS  SVs  to  the  user  will  be  on  the  horizon.  Equation 
(3.4)  shows  how  R0  is  derived,  with  a  being  the  semi-major  axis  of  the  GPS  SV  pulled 
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from  the  current  almanac  used  in  the  scenario,  and  te  the  radius  of  Earth. 


Ro  =  ^/«2-r  |  (3.4) 

Received  power  by  itself  gives  little  information  unless  compared  to  the  back¬ 
ground  noise,  i.e.  in  some  signal  to  noise  ratio.  For  GPS  applications,  C/NO  is  the 
ratio  commonly  used  because  it  refers  to  the  ratio  of  the  carrier  power  and  the  noise 
power  per  unit,  as  opposed  to  SNR  which  refers  to  the  ratio  in  a  given  bandwidth 
[9].  If  the  background  noise  temperature  is  specified,  then  the  GPS  signal  power  at 
the  receiver  can  be  converted  to  C/NO  using  Equation  (3.5). 

C/NO  =  SPrcvr  -  30  -  10  log  10{K  *  Tn)  (3.5) 

Receiver  power  is  first  converted  from  dBm  to  dBW  by  subtracting  30  dB,  K  is 
the  Boltzman’s  constant,  and  Tn  is  the  background  thermal  noise  set  at  270  Kelvin 
for  this  research.  Once  the  C/NO  is  calculated,  it  is  compared  with  a  preset  criteria 
to  determine  if  the  signal  strength  is  strong  enough  to  be  discerned  by  the  receiver. 

3.6  Criteria  for  GPS  in  View 

The  main  purpose  of  the  LOS,  Doppler,  and  C/NO  calculations  is  to  determine 
which  GPS  Pseudo-Random  Noises  (PRNs)  the  user  has  access  to.  A  connection  is 
considered  established  if  all  three  values  of  LOS,  Doppler  shift  magnitude,  and  C/NO 
meet  the  preset  criteria.  For  LOS  criteria,  the  Earth  to  Signal  range  is  required  to  be 
greater  than  the  Earth’s  radius  plus  troposphere  height.  Troposphere  is  included  into 
the  Earth’s  radius  to  avoid  the  issue  of  errors  induced  from  traveling  through  multiple 
troposphere  regions.  Doppler  shift  criteria  is  purely  receiver  based,  but  terrestrial 
receivers  have  a  height  and  velocity  restriction  in  place.  Special  firmware  can  be 
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obtained  to  remove  the  restrictions  and  even  increase  the  Doppler  search  window  up 
to  45  kHz  as  in  [7].  For  this  research,  the  Doppler  criteria  is  therefore  set  to  45  kHz. 
Finally,  the  C/NO  threshold  will  be  varied  but  kept  between  30  — 35dBHz.  This  range 
was  decided  based  on  previous  research  where  a  GPS  receiver  in  a  Geosynchronous 
Transfer  Orbit  observed  C/NO  values  ranging  from  30.41  —  42.72  dBHz  [16].  Once 
the  GPS  PRNs  meeting  all  three  thresholds  are  determined,  the  positions  and  LOS 
ranges  are  stored  off  to  be  used  in  a  position  determination  algorithm. 

3.7  Range  Noise 

Noise  is  added  to  the  perfect  LOS  ranges,  i.e.  GPS  to  user  range,  to  simulate 
a  realistic  environment.  Normally  distributed  random  numbers  are  generated  to  use 
as  errors  for  range  measurements  in  meters.  The  standard  deviation  for  the  noise  is 
derived  from  range  measurement  errors  a  GPS  user  commonly  experiences.  These 
range  measurement  errors  are  referred  to  as  UERE.  UERE  errors  are  considered 
to  be  normally  distributed  with  zero  mean.  A  standard  deviation  value  of  approxi¬ 
mately  5.36  m  was  chosen  as  the  UERE  in  this  research  and  is  described  more  in  the 
subsection  below. 

User  Range  Error. 

Range  measurement  errors,  referred  to  as  UERE,  are  the  effective  accuracy  of 
pseudorange  measurements  from  a  GPS  satellite  to  the  user.  UERE  for  a  satellite 
is  a  combination  of  each  error  source  for  that  satellite  and  is  normally  described 
as  a  one  sigma  error  in  meters.  Each  satellite’s  error  sources  are  assumed  to  be 
independent  but  the  UERE  can  be  modified  to  accommodate  correlated  error  between 
satellites  [11].  UERE  error  sources  can  be  grouped  into  two  categories:  Space  and 
User.  Space  errors  are  those  that  originate  from  the  GPS  satellite  itself.  The  space 
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error  sources  consist  of  ephemeris  error,  information  about  the  position  of  the  GPS 
SV,  and  satellite  clock  error.  User  error  sources  consist  of  troposphere,  ionosphere, 
receiver,  and  multipath  error.  Table  3.1  lists  the  current  error  values  from  the  National 
Oceanic  and  Atmospheric  Administration  as  of  April  2014  [8]. 


Table  3.1.  The  GPS  Error  Budget 


ERROR 

VALUE  (in  meters) 

Ionosphere 

4.0 

Ephemeris 

2.1 

Clock 

2.1 

Troposphere 

0.7 

Receiver 

0.5 

Multipath 

1.0 

TOTAL 

10.4 

UNCORRELATED  ERROR  (UERE) 

5.15  (rss  of  errors) 

This  research  uses  a  larger  UERE  value,  5.36  m,  than  the  value  of  5.15  m  listed  in 
Table  3.1.  Ionosphere  and  troposphere  errors  are  reduced  due  to  assumptions  made 
about  use  of  ionosphere-free  measurements  and  the  increased  radius  of  Earth.  The 
receiver  error  is  increased  to  account  for  unmodeled  error  and  the  variety  of  available 
GPS  receivers.  For  this  research,  the  UERE  value  is  not  scaled  to  account  for  variable 
path  loss  among  signals  from  multiple  SVs;  the  UERE  also  does  not  increase  due  to 
the  GPS  signal  originating  from  the  GPS  antenna  sidelobes. 

Thermal  Noise  Error. 

Additional  noise  is  added  to  each  range  measurement  based  on  the  respective 
C/NO  values  for  each  GPS  signal.  Thermal  noise  is  the  main  contributor  to  carrier 
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tracking  error  in  the  receiver’s  PLL  [11].  Carrier  tracking  error  affects  range  error 
and  is  directly  related  with  C/NO: 

c  1 

ar“  =  JTi2i 

c  is  the  speed  of  light,  fn  is  the  frequency  of  the  LI  band,  Bn  is  the  PLL  noise 
bandwidth,  Cn o  is  the  Carrier  Power  to  Noise  Density  Ratio,  and  PIT  is  the  receiver’s 
Predetection  Integration  Time  (PIT).  Before  selecting  values  for  Bn  and  PIT  (which 
are  receiver  dependent),  an  assumption  must  be  made  about  the  type  of  PLL  model 
to  be  used.  Note  this  research  bases  the  noise  calculation  off  of  the  LI  frequency  due 
to  the  assumption  that  the  LI  band  is  the  only  frequency  processed  by  the  receiver. 
This  research  effort  assumes  a  third  order  PLL  is  used  for  stability  in  a  high  dynamic 
environment.  For  a  third  order  PLL  to  be  stable,  Bn  must  be  less  than  or  equal  to 
18  Hz,  so  this  was  the  value  chosen.  PIT  values  are  normally  set  between  1  ms  and 
20  ms.  Lower  PIT  times  are  chosen  in  situations  where  high  responsiveness  is  desired, 
whereas  higher  PIT  times  are  selected  for  system  robustness  against  errors  stemming 
from  noise.  A  PIT  value  of  20  ms  is  chosen  due  to  the  user’s  six  classical  orbital 
elements  being  known  a  priori  allowing  the  orbit  to  be  predictable. 

The  ranges  with  the  noise  added  are  considered  pseudoranges  and  are  input  to 
the  position  determination  algorithm. 

3.8  Physical  Simulation 

A  second  step  to  this  research  is  data  comparison  with  position  determined  from 
a  replicated  GPS  signal.  The  signal  replication  will  be  handled  by  the  SPIRENT 
GSS800  which  produces  a  RF  signal  that  takes  into  account  propagation  loss  as  well 
as  the  GPS  transmit  antenna’s  sidelobes. 

The  SPIRENT  GSS8000  is  configured  with  the  same  scenario  information  that 
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was  input  into  the  MATLAB  GUI  to  include  start  time,  end  time,  and  relevant  GPS 
almanac.  User  position  is  generated  from  the  User  Motion  Trajectory  file  produced 
by  the  MATLAB  GUI.  Once  the  scenario  is  started,  the  SPIRENT  simulator  gener¬ 
ates  the  RF  signal  expected  to  be  received  at  the  user’s  location  in  real  time.  The 
generated  signal  is  recorded  for  approximately  36  seconds  to  ensure  one  full  frame  of 
navigation  data,  30  seconds,  is  present.  This  process  is  done  for  different  sections  in 
the  user’s  orbit.  Afterwards,  the  recorded  data  is  processed  with  an  SDR  to  calculate 
pseudoranges  for  each  timestep. 

Software  Defined  Receiver. 

The  SDR  used  for  this  research  is  also  located  within  MATLAB.  An  RF  front  end  is 
used  to  record  and  downsample  the  transmitted  signal  from  the  SPIRENT  GSS8000. 
Individual  GPS  signals  are  found  by  correlating  the  PRN  modulated  on  the  carrier 
wave  with  a  locally  generated  PRN  for  each  GPS  satellite.  Doppler  along  with  the 
carrier  wave  are  then  removed  from  each  signal,  leaving  the  navigation  data  remaining. 
GPS  transmission  time,  in  units  of  GPS  week  seconds,  can  be  determined  from  the 
navigation  message.  Transmit  time  is  then  subtracted  from  the  time  the  receiver 
first  observed  the  signal,  also  in  GPS  week  seconds.  The  difference  is  multiplied  by 
the  speed  of  light  to  generate  pseudoranges.  This  process  is  well  established  and 
more  information  can  be  found  in  [14] .  These  pseudoranges  are  then  put  through  the 
same  position  determination  algorithm  as  the  pseudoranges  from  the  supplementary 
analysis. 

3.9  Summary 

This  chapter  provides  an  overview  of  the  process  involved  in  this  research,  the 
relevant  assumptions  made,  and  the  steps  taken  to  validate  the  code  used  in  the 
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research.  Afterwards,  an  in  depth  description  of  the  supplementary  analysis  algorithm 
is  reviewed,  to  include  LOS  calculation,  Doppler  shift  calculation,  received  signal 
strength,  and  the  process  of  selecting  GPS  SVs  the  user  is  capable  of  accessing. 
Finally,  noise  generation  for  the  creation  of  pseudoranges  and  the  SDR  is  described. 
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IV.  Results  and  Analysis  of  GPS  Based  Position  Solutions 
at  Different  Orbitial  Altitudes 

4.1  Introduction 

A  total  of  seven  different  orbits  are  simulated  to  determine  the  number  of  GPS 
SVs  available  to  the  user  at  various  spots  in  its  orbit.  The  magnitude  of  errors 
expected  while  using  those  GPS  SVs  to  resolve  position  is  also  calculated.  Each 
scenario  simulation  began  on  January  26,  2015,  and  ran  for  exactly  24  hours  with  one 
minute  increments,  terminating  on  January  27th,  2015.  A  24  hour  period  is  chosen 
for  the  simulation  length,  ensuring  two  orbital  periods  of  GPS  are  present. 

The  orbits  chosen  represent  various  altitudes  above  Earth,  to  include  those  in 
LEO,  MEO,  and  one  in  GEO.  For  each  scenario,  a  breakdown  of  the  number  of 
GPS  SVs  meeting  and  failing  the  specified  criteria  is  reviewed  first.  Afterwards,  the 
DOP  and  position  error  standard  deviations  are  compared  against  the  user’s  range 
as  well  as  each  other.  Eigenvalues  of  the  errors  and  DOP  are  reviewed  to  check  for 
consistency. 

With  the  exception  of  the  GEO  scenario,  the  other  orbits  have  orbital  periods 
that  are  shorter  than  the  24  hour  scenario  duration.  This  allows  all  but  one  scenario 
to  contain  more  than  one  revolution  of  the  user’s  orbit  present  in  the  simulation. 
This  research  takes  advantage  of  the  user’s  multiple  revolutions  by  representing  some 
information  with  respect  to  the  orbital  period,  “time  since  start  of  user’s  orbital 
period”,  as  opposed  to  “time  since  start  of  scenario”.  All  figures  that  are  plotted 
using  a  time  span  less  than  1440  minutes  (24  hours)  to  be  in  reference  to  the  orbital 
period. 
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4.2  General  Parameters 


The  parameters  listed  in  Table  4.1  remain  constant  for  all  the  scenarios  presented 
in  this  research  unless  otherwise  stated. 

Table  4.1.  General  Parameters 


Parameter 

Value 

Notes 

Earth’s  Radius 

6478.1  km 

Includes  a  troposphere  height  of  100  km 

Noise  Temperature 

230  Kelvin 

Doppler  Threshold 

45  kHz 

Improved  Doppler  search  window  from  [7] 

C/N0  Threshold 

33dBHz 

Observed  values  of  30.41  —  42.91  in  GTO[16] 

GPS  Frequency 

Ll:  1575.42  MHz 

PIT 

20  ms 

PLL  Noise  Bandwidth 

18  Hz 

3rd  order  PLL 

Range  Noise  sigma 

5.357  m 

Derived  from  UERE 

Creation  of  the  different  orbital  altitudes  is  accomplished  by  varying  the  Semi- 
Major  axis,  causing  the  orbital  parameters  to  be  similar  (with  the  exception  of  the 
GEO  orbit).  As  the  orbits  are  circular,  perigee  does  not  exist  thus  the  argument  of 
perigee  is  not  defined.  Circular  orbits  also  have  no  eccentricity,  so  eccentricity  is  set 
to  0.  The  right  ascension  of  the  ascending  node  is  set  to  0°  and  the  true  anomaly 
also  starts  at  0°  for  each  scenario.  These  four  values  will  not  be  listed  in  the  tables 
to  follow. 

4.3  Low  Earth  Orbit 

Three  altitudes  were  chosen  for  the  LEO  scenario  group:  400  km,  500  km,  and 
600  km.  Each  of  the  resulting  orbits  have  a  nominal  orbital  period  of  92  min.  The 
classical  orbital  parameters  used  to  described  the  orbit  are  listed  in  Table  4.2. 


Table  4.2.  LEO  Scenarios:  Orbital  Parameters 


Altitude 

Semi-Major  Axis 

Inclination 

Period 

400  km 

6,778.137  km 

45° 

92  min 

500  km 

6,878.137  km 

45° 

94  min 

600  km 

6,978.137  km 

45° 

96  min 

Doppler  at  400km  Altitude. 

The  velocity  required  to  maintain  orbit  results  in  the  Doppler  between  the  user  and 
GPS  SVs  to  values  higher  than  that  experienced  on  Earth.  Orbital  speed  increases  as 
the  altitude  decreases.  A  user  at  400  km  altitude  will  have  Doppler  shifts  greater  than 
the  other  orbits  chosen  for  this  research  due  to  the  higher  orbital  velocity.  Figure  4.1 
shows  that  the  maximum  value  the  Doppler  can  reach  during  the  400  km  simulation 
is  around  40  kHz. 
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Figure  4.1.  Maximum  and  minimum  Doppler  by  time  at  400km 

The  Doppler  search  window  for  this  research  is  assumed  to  be  expanded  to  45  kHz 
based  on  data  found  in  [16] .  The  increased  Doppler  search  window  causes  the  velocity 
experienced  at  LEO  to  no  longer  be  a  limiting  factor  to  the  number  of  GPS  SVs  in 
view.  A  relative  speed  of  8.563  km/s  is  required  between  the  GPS  and  user  to  fail  a 
Doppler  threshold  of  45  kHz.  Failure  due  to  Doppler  shift  does  not  appear  in  any  of 
the  criteria  breakdown  figures,  such  as  Figure  4.2. 

GPS  Connection  Overview. 

Figure  4.2  shows  the  number  of  GPS  SVs  that  either  passed  or  failed  the  criteria 
for  connection  during  the  400  km  simulation. 
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Figure  4.2.  GPS  Satellite  Connection  at  400km 

For  the  entire  24  hour  (1440  min)  duration  of  the  400  km  scenario,  a  minimum  of  13 
GPS  SVs  passes  all  the  criteria  necessary  to  be  considered  usable,  as  indicated  in  blue. 
GPS  SVs  that  failed  to  meet  the  connection  criteria  did  so  because  of  obstruction 
from  Earth,  and  are  shown  in  red.  The  yellow  horizontal  line  indicates  access  to  at 
least  four  GPS  SVs  for  that  particular  time  instance.  The  times  with  access  to  less 
than  four  GPS  SVs  are  indicated  by  red  instead  of  yellow  on  the  horizontal  line. 

The  results  for  the  other  LEO  scenarios,  chosen  at  altitudes  500  km  and  600  km, 
are  similar  to  the  400  km  case. 

Figures  4.3a  and  4.3b  both  show  the  500  km  and  600  km  scenarios  having  a  mini¬ 
mum  of  13  GPS  SVs  in  view  for  the  entire  scenario  duration.  Failing  the  LOS  criteria, 
being  obstructed  by  Earth,  is  the  only  criteria  for  failure  of  any  of  the  GPS  SVs. 
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(a)  GPS  Satellite  Connection  at  500km 
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(b)  GPS  Satellite  Connection  at  600km 
Figure  4.3.  GPS  Satellite  Connection  for  LEO  Scenarios 
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DOP. 


Being  positioned  beneath  the  GPS  orbit,  along  with  the  large  number  of  GPS  SVs 
in  view,  benefited  the  user.  It  increased  the  chances  for  a  high  angle  of  separation 
existing  between  GPSs  SVs,  relative  to  the  user.  High  angle  of  separation  directly 
contributes  to  a  lower  DOP.  DOP  is  calculated  based  on  the  true  position  of  the 
user  and  the  positions  of  the  GPS  SVs  meeting  all  three  access  criteria.  GDOP 
experienced  by  the  user  in  the  400  km  altitude  scenario  remained  under  two  for  the 
entire  duration.  GDOP  and  PDOP  for  the  24  hour  scenario  are  shown  in  Figure  4.4 
as  a  function  the  user’s  orbital  period. 
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Figure  4.4.  Comparison  of  the  GDOP  and  PDOP  vs  time  at  400km 

In  Figure  4.4,  GDOP  is  greater  in  magnitude  but  follows  the  same  trend  as  PDOP. 
Due  to  the  similarities,  GDOP  is  the  DOP  metric  this  research  uses  to  since  the  others 
will  only  have  improved  values  relative  to  GDOP.  The  red  line  in  Figure  4.4  is  a 
visual  representation  of  a  DOP  threshold  of  20.  Any  DOP  value  above  the  threshold 
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is  considered  to  be  poor  and  unusable  for  a  terrestrial  user  based  on  [15].  The  GDOP 
values  for  this  scenario  are  well  below  this  threshold. 

Similar  to  the  400  km  scenario,  both  500  km  and  600  km  GDOP  values  also  remain 
well  under  two  for  the  entirety  of  the  simulation  duration.  GDOP  for  all  three  LEO 
altitudes  is  shown  in  Figure  4.5.  This  shows  that  in  addition  to  no  significant  loss 
in  the  number  of  SVs  available,  there  is  no  significant  change  in  geometry  over  this 
altitude  range. 


Figure  4.5.  GDOP  vs  Time  at  LEO 


Position  Errors. 

Position  errors  presented  in  this  chapter  refer  to  the  distances  between  true  user 
positions  and  the  results  of  the  Newton-Raphson  least  squares  calculations  over  time. 
Errors  in  position  estimates  for  the  user  at  400  km  altitude  are  similar  to  those  experi¬ 
enced  by  terrestrial  GPS  users.  The  overall  errors  are  below  10  meters  for  a  majority 
of  the  user’s  orbit,  with  the  exception  of  a  few  instances  where  the  error  rises  above 
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10  meters.  Figure  4.6  shows  a  trace  of  the  user’s  orbit  for  one  Monte  Carlo  instance 
for  the  entire  scenario  duration.  The  color  of  the  dots  represent  the  magnitude  of  the 
position  estimation  error,  true  user  position  compared  against  the  position  calculated 
from  the  Newton-Raphson  algorithm.  The  percentages  represented  in  the  legend  are 
calculated  based  on  the  current  Monte  Carlo  run  represented. 

Error  Mag  <  1 0m  :  98.5% 

1 0m  >=  Error  Mag  <  1 00m  :  1 .5% 


Figure  4.7  shows  the  standard  deviation  of  the  overall  position  error  by  orbital 
period  for  all  of  the  LEO  scenarios.  For  the  duration  of  the  user’s  orbital  period,  the 
rms  error  standard  deviation  remains  under  3  meters  for  all  three  LEO  scenarios. 
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Figure  4.7.  RMS  Error  STD  for  LEO  scenarios 


A  breakdown  of  the  position  error  standard  deviation  for  each  axis  in  the  400  km 
scenario  is  shown  in  Figure  4.8.  Over  1000  Monte  Carlo  trials,  the  error  standard 
deviation  for  all  three  axes  each  remain  under  six  meters. 


46 


Figure  4.8.  Error  STD  by  axis  vs  User  Range  at  400km 


Similar  to  the  400  km  scenario,  position  errors  for  the  500  km  and  600  km  alti¬ 
tude  scenarios  also  remain  under  10  m  with  standard  deviations  remaining  under  3  m. 
Figures  4.9a  and  4.9b  display  the  position  error  over  the  orbit. 
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Figure  4.9.  Position  Error  Over  Orbital  Trajectory  for  LEO  Scenarios 
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Applicability  of  DOP. 


Another  goal  of  this  research  is  to  determine  if  DOP  is  applicable  for  predicting 
precision  for  non-terrestrial  users  of  GPS.  For  the  LEO  scenarios,  the  magnitude 
of  both  DOP  and  position  errors  are  consistent  with  those  of  a  terrestrial  user.  To 
determine  consistency  between  DOP  and  position  error,  it  is  necessary  to  compare 
two  results  against  each  other.  In  order  to  do  so  DOP  had  to  be  weighted  with 
the  range  errors  expected  for  each  GPS  SV  in  view.  Weighting  was  accomplished 
by  multiplying  the  UERE  with  DOP  to  obtain  predicted  position  sensitivity.  The 
comparison  of  the  400  km  scenario’s  DOP  and  position  error  are  shown  in  Figure 
4.10. 
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Figure  4.10.  Position  Error  STD  by  axis  vs  DOP  at  400km 


Plotting  the  position  error  standard  deviation  against  each  axis’s  DOP  helps  de¬ 
termine  if  the  linear  relationship  between  position  error,  DOP,  and  range  error  exists. 
Figure  4.11a  shows  the  standard  deviation  of  position  error  plotted  again  the  respec- 
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tive  DOP  values  for  each  axis.  Figure  4.11b  shows  the  DOP  vs  position  error  STD 
broken  up  into  the  three  axes. 
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(a)  DOP  vs  Combined  Position  Error  STD  at  400km 
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(b)  DOP  vs  Position  Error  STD  for  each  Axis  at  400km 
Figure  4.11.  DOP  vs  Position  Error  STD  at  400km 


50 


A  table  of  the  linear  fit  values  for  all  scenarios  can  be  found  after  the  GEO  section 
in  Table  4.5.  Applying  a  linear  fit  to  the  combined  data  points  produces  a  line  with 
a  slope  of  .00535  km  and  a  offset  of  approximately  0.  Each  axis  also  shares  a  similar 
STD  to  DOP  linear  fit  equation,  where  the  slope  is  around  .0053  km.  This  value 
is  consistent  with  the  previously  set  range  noise  sigma  of  5.357m  derived  from  the 
UERE,  proving  that  the  linear  relationship  still  exists. 

The  linear  relationship  also  holds  true  for  the  remaining  LEO  simulations  as  seen 
in  Figures  4.12a  and  4.12b. 
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Figure  4.12.  DOP  vs  Position  Error  STD  for  LEO  Scenarios 
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To  see  if  the  DOP  and  position  errors  are  aligned  within  the  same  space,  the 
Eigenvalues  for  position  error  and  DOP  are  compared.  Figure  4.13  compares  the 
Eigenvalues  between  GDOP  and  the  error  RMS  standard  deviation,  while  Figure 
4.14  shows  the  comparison  for  each  axis  over  the  scenario  duration. 
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Figure  4.13.  Error  vs  DOP  Eigenvalue  Comparison  at  400km 
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Figure  4.14.  Eigenvalue  Comparison  by  axis  at  400km 


Both  plots  shows  that  the  position  errors  still  follow  the  respective  axis’s  DOP.  Sim¬ 
ilar  results  are  seen  with  the  500  km  and  600  km  scenarios  in  Figures  4.15  and  4.16. 
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Figure  4.15.  Error  vs  DOP  Eigenvalue  Comparison  at  500km 
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Figure  4.16.  Error  vs  DOP  Eigenvalue  Comparison  at  600km 


Based  on  the  result  of  these  three  LEO  altitudes,  DOP  is  still  a  valid  metric  for 
predicting  position  error  sensitivity  based  on  GPS  SV  geometry  for  a  user  in  LEO. 
In  addition,  the  software  is  validated  to  perform  the  expected  calculations  according 
to  the  sensitivity  analysis/measurement  theory. 
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4.4  Medium  Earth  Orbit 


Three  altitudes  were  chosen  for  the  MEO  scenario  group:  10100  km,  15150  km, 
and  25000  km.  These  altitudes  were  chosen  relative  to  the  GPS  constellations  alti¬ 
tude.  The  orbits  are  at  \  GPS  altitude,  |  GPS  altitude,  and  slightly  above  the  GPS 
constellation  altitude.  The  classical  orbital  parameters  for  these  orbits  are  listed  in 
Table  4.3. 


Table  4.3.  MEO  Scenarios:  Orbital  Parameters 


Altitude 

Semi-Major  Axis 

Inclination 

Period 

10, 100  km 

16,478.137  km 

45° 

351  min 

15,150  km 

21,528.137km 

45° 

524  min 

25,000  km 

31,378.137km 

45° 

922  min 

Comparing  the  number  of  GPS  SVs  in  view  for  the  LEO  and  10100  km  scenario 
brings  about  an  unexpected  result.  Figure  4.17  indicates  that  a  receiver  at  10100  km 
will  have  double  the  available  GPS  SVs  than  those  available  at  LEO. 
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Figure  4.17.  GPS  Satellite  Connection  at  10100km 


Earth  becomes  less  of  an  obstruction  the  further  away  the  user  is.  However,  this 
also  increases  the  likelihood  of  being  outside  of  the  GPS  transmit  antenna’s  main 
beam  area.  The  large  number  of  GPS  SVs  available  to  a  user  at  an  orbital  altitude  of 
10100  km  is  due  the  effect  of  GPS  antenna’s  side  lobes.  Figure  4.18  shows  the  result 
of  running  the  same  scenario  again  but  only  considering  effects  from  the  main  beam 
of  the  GPS  antenna. 
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Figure  4.18.  Satellites  in  View  only  considering  GPS’s  main  beam  at  10,100km 


In  this  variation  of  the  scenario,  the  user  was  not  in  the  main  beam  path  of  four 
or  more  GPS  SVs  at  any  point  in  the  24  hour  simulation  duration.  The  remainder  of 
this  research  will  only  consider  the  cases  that  use  both  the  GPS  antenna  main  beam 
and  side  lobes. 

The  number  of  GPS  SVs  in  view  drops  as  the  orbital  altitude  of  the  user  increases 
due  to  lower  signal  strength.  This  is  seen  when  looking  at  Figures  4.19  and  4.20. 
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Figure  4.19.  GPS  Satellite  Connection  at  15,150km 


Unlike  the  LEO  altitude  scenarios,  the  majority  of  connection  failures  is  not  from 
Earth  obstruction.  Rather  the  predominate  criteria  failures  originate  from  low  signal 
strength,  when  C/N0  is  less  than  the  receiver  threshold.  This  is  more  pronounced  in 
Figure  4.20. 
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Figure  4.20.  GPS  Satellite  Connection  at  25,000km 


There  is  even  a  point  in  time  where  the  user  at  25000  km  altitude  is  “in  view”  of 
only  four  GPS  SVs.  Even  though  more  SVs  have  line  of  sight  to  the  user,  the  distances 
are  greater.  The  increased  distance  causes  the  signals  originating  from  GPS  side  lobes, 
which  started  with  low  gain,  to  fall  under  the  receiver’s  C/N0  threshold  at  time  of 
reception. 
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DOP. 


DOP  for  the  10100  km  and  15150  km  altitude  scenarios  are  both  around  one. 
Different  results  are  seen  for  the  25000  km  orbital  altitude  scenario.  Figure  4.21 
displays  all  three  scenario’s  GDOP  for  comparison. 


Figure  4.21.  GDOP  vs  Time  at  MEO 


Unlike  the  10100  km  and  15150  km  altitudes,  the  25000  km  scenario  is  above  the 
GPS  constellation.  The  location  of  available  signals  is  changed  when  the  user’s  alti¬ 
tude  surpasses  that  of  the  GPS  constellation,  thus  affecting  the  GDOP.  The  values 
of  the  25000  km  scenario’s  GDOP  generally  stay  around  two  but  does  occasionally 
rise  past  10. 
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Position  Errors. 


The  effect  the  number  and  location  of  available  GPS  SVs  has  on  position  can 
be  viewed  when  comparing  the  position  errors  for  altitudes  10100  km,  15150  km,  and 
25000  km.  Figures  4.22,  4.23,  and  4.24  present  one  realization  of  the  position  errors 
expected  for  the  user’s  orbit  for  each  MEO  scenaro. 

Error  Mag  <  1 0m  :  99.9% 

10m  >=  Error  Mag  <  100m  :  0.1% 

x  104 
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Z(km)  ^  Z(km) 


Error  Mag  <  10m  :  99.7% 

10m  >=  Error  Mag  <  100m  :  0.3% 


Error  Mag  <  1 0m  :  46.4% 

10m  >=  Error  Mag  <  100m  :  53.5% 
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A  majority  of  the  position  errors  experienced  at  an  altitude  above  the  GPS  con¬ 
stellation  were  between  10  m  and  100  m,  as  opposed  to  under  10  m  for  the  altitudes 
below  the  constellation.  At  altitudes  above  the  GPS  constellation,  the  user  has  to 
rely  more  on  GPS  signals  originating  from  the  opposite  side  of  Earth.  This  not  only 
decreases  signal  strength,  it  also  decreases  the  relative  angle  between  GPS  SVs  avail¬ 
able  to  the  user.  The  reliance  on  these  signals  with  increasingly  singular  direction  of 
origin  causes  the  positional  error  to  also  be  along  that  same  direction  as  shown  in 
Figure  4.25.  The  magnitudes  of  DOP  and  the  position  errors  have  been  amplified  by 
a  factor  of  4  *  104  to  provide  a  visual  representation  along  the  orbit’s  path. 


Figure  4.25.  3D  Visualization  of  DOP  and  Error  at  25,000km 
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The  position  error  indicated  by  the  red  dots  in  Figure  4.25  are  more  elongated 
compared  to  those  of  the  10100  km  case  in  Figure  4.26. 


Figure  4.26.  3D  Visualization  of  DOP  and  Error  overlap  at  10,100km 


Figure  4.27  displays  a  3<r  representation  of  both  10100  km  and  25000  km’s  weighted 
DOP  at  the  very  beginning  of  the  simulation.  Both  positions  are  aligned  over  the 
same  location  on  Earth  at  this  time  instance. 
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Figure  4.27.  3cr  DOP  Comparison  at  Scenario  Start 


Applicability  of  DOP. 

The  MEO  scenarios  above  and  below  the  GPS  constellation  exhibited  the  linear 
relationship  between  DOP  and  position  error  similar  to  the  LEO  scenarios.  Altitudes 
10100  km  and  25000  km  are  shown  in  Figures  4.28  and  4.29  respectively. 


66 


Position  Error  STD  (km) 


xIO'3 


Figure  4.28.  DOP  vs  Position  Error  STD  at  10,100km 


Figure  4.29.  DOP  vs  Position  Error  STD  at  25,000km 
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4.5  Geostationary  Orbit 


A  Geostationary  orbit  with  an  altitude  of  35788  km  is  chosen  as  the  last  scenario 
to  test.  The  orbital  parameters  are  listed  in  Table  4.4 

Table  4.4.  GEO  Scenarios:  Orbital  Parameters 


Altitude 

Semi-Major  Axis 

Inclination 

RAAN 

35,788  km 

16,478.137  km 

45° 

1436  min 

At  this  distance  obtaining  an  abundance  of  GPS  SVs  “in  view”  is  difficult  due 
to  the  distance  from  the  GPS  constellation.  Similar  to  the  25000  km  scenario,  the 
limiting  factor  is  low  signal  strength  as  seen  in  Figure  4.30.  Since  users  at  higher  orbits 
depend  on  the  low  gain  GPS  antenna’s  side  lobes,  the  increased  distance  causes  the 
signals  to  fall  below  the  receiver’s  signal  strength  threshold.  Unlike  the  LEO  and  MEO 
altitudes,  the  number  of  GPS  SVs  in  view  for  the  GEO  scenario  never  reaches  10. 
There  are  multiple  instances  where  the  minimum  threshold,  four  GPS  SVs  present, 
is  not  met  for  position  and  DOP  calculation. 
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Figure  4.30.  GPS  Satellite  Connection  at  GEO 


DOP. 

The  few  number  of  GPS  SVs  in  view  all  originate  from  a  similar  general  direction, 
causing  the  DOP  to  worsen.  At  no  point  during  the  scenario  duration  does  the  GDOP 
fall  under  the  threshold  of  20,  shown  in  Figure  4.31. 
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Figure  4.31.  GDOP  vs  Time  at  GEO 

If  these  GDOP  values  were  received  on  Earth,  they  would  be  considered  poor  and 
possibly  advised  to  be  discarded  [15]. 

Position  Error. 

Unlike  the  LEO  and  GEO  scenarios,  the  position  error  for  the  GEO  altitude  is 
rarely  below  10  m.  When  the  GEO  position  is  able  to  be  calculated,  the  position  error 
ranges  anywhere  from  under  10  m  to  over  1km  as  shown  in  Figure  4.32. 
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Figure  4.32.  3D  Visualization  of  RMS  Error  at  GEO 


71 


Applicability  of  DOP. 


The  linear  relationship  between  DOP  and  the  position  error  is  still  present  for 
the  GEO  altitude  even  though  the  magnitude  of  the  position  errors  have  increased. 
Figure  4.33  displays  this  relationship.  The  linear  fit  values  for  each  scenario  can  be 
found  in  Table  4.5. 


Figure  4.33.  DOP  vs  Position  Error  STD  at  GEO 
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Table  4.5.  Linear  Fit  Model  for  Error  to  DOP  Relationship.  Linear  equations  takes  the 
form:  Position  Error  STD  =  Linear  Fit  *  DOP  +  Offset.  The  linear  fit  factor  in  the 
figures  above  are  in  units  of  km  while  the  table  converts  it  to  units  of  m  for  comparison 
with  the  UERE 


Scenario 

Linear  Fit 

Offset 

||  UERE  -  Fit  || 

R2 

400  km 

5.349 

3.573e-06 

0.0082 

.999 

500  km 

5.358 

-1.589e-06 

0.0007 

.999 

600  km 

5.375 

-8.274e-06 

0.0177 

.999 

10100  km 

5.369 

-3.709e-06 

0.0117 

.999 

15150  km 

5.371 

-6.733e-06 

0.0138 

.999 

25000  km 

5.359 

-2.168e-07 

0.0018 

.999 

GEO 

5.256 

-0.0029 

0.1012 

.999 

When  overlaying  the  position  errors  and  DOP  values  onto  the  user’s  orbit,  it  can  be 
shown  that  the  majority  of  the  error  is  in  the  direction  from  which  the  GPS  signals 
arrive.  The  errors  in  Figure  4.34  are  more  pronounced  compared  to  the  25000  km 
scenario  whose  altitude  is  also  above  the  GPS  constellation.  In  Figure  4.34,  the 
scaling  factor  applied  for  visualization  purposes  is  |  of  the  value  used  for  the  MEO 
scenarios  in  Figures  4.25  and  4.26. 
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DOP  and  Error  vs  Orbit  Position  (ECI) 


4.6  Visual  Summary  of  GPS  Connections 

At  the  beginning  of  the  LEO  and  MEO  scenarios,  each  user’s  relative  position 
over  Earth  is  the  same.  The  distance  away  from  Earth  is  represented  as  different 
magnitudes  along  the  ECI  X  axis.  At  the  scenario  start  time,  each  GPS  SV  with  a 
valid  connection  to  the  user  is  calculated  and  plotted  for  scenarios  500  km,  15150  km, 
and  25000  km.  Figures  4.35  -  4.37  provide  a  visualization  of  how  the  connected  GPS 
geometry  changes  with  altitude.  The  GEO  simulation  does  not  start  at  the  same 
relative  position  over  Earth.  However  Figure  4.38  provides  an  example  when  four 
GPS  SVs  are  in  view. 
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Figure  4.35.  GPS  in  View  at  Scenario  Start  at  500km 
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Figure  4.36.  GPS  in  View  at  Scenario  Start  at  15,150km 
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Figure  4.37.  GPS  in  View  at  Scenario  Start  at  25,000km 
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Figure  4.38.  GPS  in  View  at  Scenario  Start  at  GEO 
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4.7  RF  Signal  Generation 


A  GPS  RF  signal  generator  was  operated  in  an  attempt  to  validate  the  results  of 
the  MATLAB  supplementary  analysis.  The  four  scenarios  chosen  to  be  loaded  into 
the  SPIRENT  GSS8000  signal  generator  were  the  400  km,  15150  km,  25000  km,  and 
GEO  scenarios.  The  resulting  RF  signal  for  each  scenario  was  recorded  for  37  second 
and  then  processed  by  a  SDR. 

GPS  SV  Acquisition. 

After  being  recorded,  the  RF  signal  is  processed  to  determine  the  number  of  GPS 
signals  present.  Simplistically  speaking,  the  recorded  signal  is  correlated  against  mul¬ 
tiple  versions  of  the  reference  signal  (accounting  for  shifts  in  Doppler  and  PRN  code). 
This  correlation  process  is  repeated  for  a  preset  number  of  PRN  code  iterations  (civil¬ 
ian  signal’s  PRN  repeats  every  millisecond).  The  multiple  versions  of  the  reference 
signal  that  were  correlated  are  converted  to  power  and  then  analyzed  by  compar¬ 
ing  the  power  against  one  another.  A  signal  is  considered  present  if  the  power  of 
the  strongest  correlation  compared  to  that  of  the  2nd  strongest  is  above  a  certain 
threshold. 

The  ratio  between  the  two  strongest  signals,  known  as  the  acquisition  ratio,  is 
a  scale  to  determine  the  likelihood  of  a  specific  GPS  SV  being  present.  A  GPS  SV 
is  considered  present  if  its  acquisition  ratio  is  above  the  preset  threshold  of  3  for 
this  research.  An  acquisition  threshold  of  three  requires  the  power  of  the  highest 
correlation  to  be  a  minimum  of  three  times  greater  than  the  next  highest  correlation. 
Figures  4.39  to  4.42  are  the  results  of  processing  the  recorded  RF  signal  from  the  four 
scenarios. 
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Figure  4.39.  GPS  Signal  Acquisition  Strength  at  400km 


PRN 

Figure  4.40.  GPS  Signal  Acquisition  Strength  at  15,150km 
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At  400  km,  Figure  4.39,  eight  GPS  SVs  are  present  and  each  have  an  Acquisition 
ratio  of  15  or  greater.  The  15, 150  km  scenario  in  Figure  4.40  has  13  GPS  SVs  present 
but  only  three  are  above  an  Acquisition  ration  of  15.  Five  of  the  signals  are  barely 
over  the  threshold,  meaning  the  signal  strength  is  so  low  as  to  be  almost  be  indis¬ 
tinguishable  from  noise.  Picking  up  signals  that  originate  from  the  side  lobes  of  the 
GPS  antenna  would  produce  these  types  of  results. 
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Figure  4.41.  GPS  Signal  Acquisition  Strength  at  25,000km 
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Figure  4.42.  GPS  Signal  Acquisition  Strength  at  GEO 


The  processed  RF  data  is  showing  14  GPS  SVs  in  view  for  the  25,  000  km  scenario 
and  15  for  the  GEO  scenario,  Figures  4.41  and  4.42  respectively.  The  signals  in  the 
25,000  km  RF  data  have  stronger  Acquisition  thresholds  than  the  15,150  km  data. 
Results  from  the  GEO  scenario  was  vastly  different  than  the  MATLAB  analysis, 
showing  well  over  the  predicted  maximum  of  10  SVs  in  view.  Having  four  PRNs  that 
distinctly  standout  compared  to  the  rest  was  predicted  for  the  GEO  case,  but  the 
sheer  number  of  GPS  signals  found  was  unexpected  so  there  are  some  reservations 
about  the  RF  results. 

One  method  to  validate  the  RF  signal  would  be  to  run  the  same  scenario  multiple 
times  and  record  the  resulting  position  errors.  Afterwards,  analyze  the  position  errors 
for  a  Gaussian  trend  that  matches  the  expected  error  ellipse.  The  results  of  the  RF 
signal  can  also  be  checked  by  analyzing  the  positions  of  the  GPS  SVs,  whose  PRN 
was  observed  in  the  RF  signal,  with  the  location  of  the  receiver  to  ensure  that  signal 
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reception  is  possible. 


Low  Acquisition  Ratio. 

Choosing  an  appropriate  acquisition  threshold  is  important  for  limiting  the  amount 
of  weak  signals.  If  a  signal  is  too  weak,  then  it  may  be  difficult  to  read  the  navigation 
data  in  the  signal.  Errors  could  also  arise  due  to  the  noise’s  effect  on  the  signal  as 
well.  To  illustrate  the  difference  in  strength,  PRNs  10  and  18  were  from  the  GEO  RF 
data  were  processed  and  the  tracking  results  are  shown  with  Figures  4.43  to  4.46. 

Plotting  the  two  outputs  of  the  RF  signal  after  it’s  been  correlated  and  accumu¬ 
lated  gives  information  on  how  well  the  signal  is  being  tracked.  If  the  receiver  is  not 
tracking  the  GPS  signal  or  no  signal  is  present,  then  all  the  data  points  will  be  at  the 
center  of  the  plot.  If  the  signal  is  strong  and  the  receiver  is  tracking  both  the  GPS 
signal’s  frequency  and  phase,  then  the  data  points  will  be  concentrated  in  two  small 
circles  far  away  from  the  center  as  shown  in  Figure  4.43  for  PRN  10. 
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Figure  4.43.  GEO  Signal  Tracking  for  PRN  10 


Even  while  able  to  be  tracked,  PRN  18’s  signal  strength  is  not  as  a  strong  so  the 
data  points  will  be  closer  to  the  center  of  the  plot,  as  shown  in  Figure  4.44. 
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Prompt  Output 
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Figure  4.44.  GEO  Signal  Tracking  for  PRN  18 

If  the  same  data  is  plotted  against  time,  data  bit  transitions  can  be  observed  if 
the  signal  strength  is  strong  enough.  PRN  10  is  plotted  in  Figure  4.45  as  an  example 
of  a  strong  signal  case  while  PRN  18  is  plotted  in  Figure  4.46  is  an  example  of  the 
contrary. 
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Early  13  and  G3 


Figure  4.45.  Tracked  and  Correlated  Signal  for  PRN  10  at  GEO 
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Figure  4.46. 


Tracked  and  Correlated  Signal  for  PRN  18  at  GEO 
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4.8  Summary 


This  chapter  compared  the  position  errors  and  resulting  DOP  for  orbital  altitudes 
that  positioned  the  user  from  LEO  to  GEO.  It  is  shown  that  the  altitudes  lower  than 
the  GPS  constellation  experienced  position  errors  under  10  m  when  utilizing  the  GPS 
antenna  sidelobes.  As  the  orbital  altitude  increases  from  LEO  until  about  |  of  the 
GPS  constellation,  the  number  of  GPS  SVs  in  view  also  increased.  Increasing  the 
altitude  from  that  location  causes  the  number  of  available  GPS  to  decrease  due  to 
lower  signal  strength  received.  Access  to  more  than  four  GPS  SVs  while  at  GEO 
became  difficult  and  the  errors  greatly  increased.  The  errors  accrued  are  in  the  same 
direction  as  the  incoming  GPS  signals.  The  linear  relationship  between  DOP  and 
position  errors  existed  for  all  the  scenarios  simulated  in  the  research.  This  shows  that 
DOP  is  still  a  valid  metric  for  predicting  precision  levels  in  altitudes  exceeding  LEO. 

Results  from  the  RF  signal  processed  by  the  SDR  showed  trends  similar  to  the 
results  of  the  MATLAB  analysis.  The  signals  uncovered  in  the  400  km  LEO  case 
are  all  distinctly  above  the  noise  floor  (acquisition  level  of  the  non-present  signals) 
and  acquisition  threshold.  At  15150  km  altitude,  the  processed  RF  contained  more 
GPS  SVs  than  the  LEO  signal,  though  a  portion  of  these  signals  are  close  to  the 
noise  floor  and  barely  passed  the  acquisition  threshold.  The  scenarios  above  the 
GPS  constellation  showed  an  unexpected  number  of  signals  present.  Based  on  this 
discrepancy,  geometric  analysis  may  be  insufficient  for  analysis,  especially  at  GEO. 
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V.  Conclusion 


5.1  Contributions 

The  research  outlined  in  this  thesis  involved  developing  a  process  to  evaluate  ex¬ 
pected  GPS  accuracy  tailored  to  the  concerns  of  satellite  orbits.  The  process  described 
in  Chapter  III  involved  modeling  the  orbits  of  the  user  and  GPS  constellation  over  a 
24  hour  period  to  obtain  position  and  velocity  information.  This  information  is  pro¬ 
cessed  and  compared  against  preset  criteria  thresholds  on  line  of  sight,  Doppler  shift 
magnitude,  and  received  signal  strength  to  determine  which  GPS  SVs  were  accessible 
to  the  user.  Pseudoranges  are  created  for  the  accessible  GPS  SVs  by  calculating  the 
true  GPS  to  user  range,  then  adding  noise  based  on  User  Equivalent  Range  Error 
and  the  C/NO  for  each  signal.  These  pseudoranges  are  used  to  evaluate  position  esti¬ 
mation  accuracy  and  the  applicability  of  Dilution  of  Precision  for  users  in  the  space 
environment. 

Seven  user  orbits  ranging  from  LEO  up  to  GEO  are  simulated,  with  two  of  the 
orbits  having  an  altitude  above  the  GPS  constellation.  All  the  scenarios  below  the 
GPS  constellation  maintained  access  to  a  minimum  of  13  GPS  SVs  throughout  for 
the  entire  24  hour  simulation  duration.  The  number  of  GPS  PRNs  steadily  dropped 
as  the  altitude  increased  for  the  scenarios  above  the  GPS  constellation.  The  GEO 
scenario  encountered  multiple  instances  where  3  or  less  GPS  PRNs  are  available. 

The  metrics  discussed  in  this  section  are  based  off  of  a  total  of  1,000  Monte  Carlo 
executions,  with  1,441  time  steps/data  points  for  each  Monte  Carlo  run.  Each  time 
instance  is  treated  as  a  trial,  resulting  in  1, 441,  000  samples  of  position  rms  error  over 
the  Monte  Carlo  runs  for  each  scenario.  Error  percentages  are  calculated  by  finding 
the  number  of  trials  out  of  1,441,000  data  points  that  meet  the  desired  criteria,  e.g. 


under  10  m. 


The  results  of  1,000  Monte  Carlo  simulation  runs  for  position  estimation  are  as 
follows:  99%  the  rms  position  errors  stayed  under  10  m  for  the  LEO  scenarios.  The 
position  error  for  both  MEO  scenarios  below  the  GPS  constellation’s  altitude  rose 
above  10  m  for  0.2%  of  the  total  trails.  These  results  are  specific  to  the  scenarios  where 
the  GPS  antenna’s  side  lobes  are  modeled.  Running  a  MEO  scenario  without  the 
side  lobe  model  resulted  in  less  than  four  GPS  SVs  available  for  the  entire  simulation 
duration.  As  expected,  the  two  orbits  above  the  GPS  constellation  showed  the  worst 
results  out  of  the  group.  In  the  MEO  scenario  with  altitude  of  25,  000  km,  46%  of 
position  errors  fell  under  10  m  and  53.5%  are  between  10  m  and  100m.  The  GEO 
scenario  has  the  worst  results:  31.2%  are  unable  to  calculate  a  solution,  0.4%  of  the 
errors  are  below  10  m,  23%  are  between  10  m  and  100  m,  and  45.4%  of  the  errors  are 
greater  than  100  m. 

DOP  is  calculated  for  every  scenario  and  compared  against  position  errors.  The 
linear  relationship  between  DOP,  measurement  range  error,  and  position  error  is 
shown  to  exist  for  all  scenarios.  The  eigenvalue  trends  of  both  DOP  and  position 
error  standard  deviation  are  consistent  with  one  another.  DOP  is  still  applicable  for 
use  when  calculating  position  with  space-borne  GPS  receivers. 

5.2  Recommendations  for  Future  Research 

There  are  a  few  recommendations  for  future  research  on  this  topic.  Each  of  these 
recommendations  revolves  around  adding  more  fidelity  to  the  model: 

•  Fully  implement  a  SDR  to  not  only  track  the  GPS  signals,  but  also  provide 
position  estimates  by  decoding  the  RF  signals  from  GPS  signal  simulator. 

•  Comparing  the  results  of  different  position  estimation  algorithms  currently  used 
in  navigation  solutions,  such  as  a  Kalman  Filter,  against  those  resulting  from 
the  Newton- Raphson  method. 
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•  Consideration  of  non-isotropic  receiver  antenna  models  and  the  effect  on  existing 
errors  due  to  antenna  placement  and  orientation  on  a  receiver  satellite. 

•  Consideration  of  multi-GNSS  receivers  and  regional  space-based  augmentation 
systems. 

5.3  Summary 

This  chapter  summarizes  the  results  of  the  thesis  research  effort  as  well  as  pre¬ 
sented  future  work  opportunities.  Chapter  1  provides  the  motivation  for  researching 
on-board  satellite  position  estimation  based  on  GPS  measurements.  Chapter  2  gave 
insight  to  the  characteristics  of  the  potential  orbits  for  the  GPS  receiver  to  reside 
and  current  GNSS  constellations.  Errors  affecting  the  GNSS  signal  and  known  tech¬ 
niques  for  estimating  position  from  those  noisy  signals  was  also  presented  in  Chapter 
2.  Chapter  3  presents  the  methodology  behind  creating  a  model  that  incorporates 
the  location  and  velocities  of  GPS  SVs  and  the  user,  errors  affecting  each,  and  trans¬ 
mitter/receiver  models  to  simulate  expected  position  error.  Chapter  4  shows  that 
the  GPS  and  the  Newton  least  squares  algorithm  for  position  estimation  will  result 
in  errors  under  10  m  for  altitudes  under  that  of  the  GPS  constellation.  For  altitudes 
above  the  GPS  constellation,  the  Newton  least  squares  method  no  longer  results  in 
errors  under  10  m.  It  also  shows  that  the  relationship  for  DOP  and  position  error 
still  holds  for  all  these  altitudes. 
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Appendix  A.  Additional  Figures 
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Figure  1.1.  DOP  vs  Position  Error  STD  at  15,150km 
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GPS  in  View  -  25000km 
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Figure  1.3.  GPS  in  View  Showing  Side  Lobe  Connection 
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List  of  Acronyms 


C/NO 

Carrier  Power  to  Noise  Density  Ratio 

dB 

Decibels 

DOP 

Dilution  of  Precision 

ECEF 

Earth  Centered  Earth  Fixed 

ECI 

Earth  Centered  Inertial 

ENU 

East-North-Up 

GDOP 

Geometric  Dilution  of  Precision 

GEO 

Geostationary  Orbit 

GLONASS  Global’naya  Navigatsionnaya  Sputnikovaya  Sistema 


GNSS 

Global  Navigation  Satellite  System 

GPS 

Global  Positioning  System 

GTO 

Geostationary  Transfer  Orbit 

GUI 

Graphic  User  Interface 

HEO 

Highly  Elliptical  Orbit 

ICD 

Interface  Control  Document 

KF 

Kalman  Filter 

LEO 

Low  Earth  Orbit 

LOS 

Line  of  Sight 
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MEO 

Medium  Earth  Orbit 

PDOP 

Position  Dilution  of  Precision 

PIT 

Predetection  Integration  Time 

PLL 

Phase  Lock  Loop 

PNT 

Position,  Navigation,  and  Timing 

PRN 

Pseudo- Random  Noise 

RF 

Radio  Frequency 

SDR 

Software  Defined  Radio 

SNR 

Signal  to  Noise  Ratio 

STK 

Systems  Tool  Kit 

SV 

Space  Vehicle 

TLE 

Two-Line  Element 

UERE 

User  Equivalent  Range  Error 
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